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A well  known  test  for  uniformity  on  the  circle  given  by  B.  Ajne  has 
been  extended  to  the  hypersphere  by  numerous  authors.  However,  under 
the  hypothesis  of  uniformity  the  asymptotic  distribution  of  Ajne's  test 
statistic  on  the  hypersphere  is  somewhat  difficult  to  work  with  since  it 
involves  an  infinite  series  of  independent,  weighted,  central,  chi- 
square  random  variables.  Presented  here  is  a new  test  statistic  for 
testing  uniformity  on  the  hypersphere  which  has  a normal  asymptotic 
distribution  under  both  the  null  hypothesis  of  uniformity  and 
rotationally  symmetric,  unimodal  alternatives. 

Relations  of  the  new  test  statistic  to  Ajne's  generalized  test 
statistic  and  Rayleigh's  test  statistic  are  given.  The  new  test  is 
shown  to  have  certain  optimal  properties  against  Langevin,  generalized 
cardioid  and  other  selected  alternatives. 

In  order  to  discuss  rotationally  symmetric,  unimodal  densities  on 
the  hypersphere,  representations  in  terms  of  Gegenbauer  polynomials  are 
considered.  The  coefficients  for  these  representations,  called 

vi 


Gegenbauer  transforms,  are  given  for  many  well  known  densities.  Also, 
for  densities  represented  by  means  of  Gegenbauer  transforms,  a new 
graphical  curve  fitting  aid  is  presented. 

In  the  course  of  investigating  Gegenbauer  polynomials,  a 
standardization  proposed  by  J.  G.  Saw  for  the  Gegenbauer  polynomials  and 
a weight  function  differing  by  a constant  factor  from  the  usual  weight 
function  for  Gegenbauer  polynomials  are  shown  to  simplify  certain 
results.  tony  formulas  for  Gegenbauer  polynomials  are  collected  and 
written  in  terms  of  Saw's  standardization  and  this  modified  weight 


function 


CHAPTER  1 

INTRODUCTION  AND  PRELIMINARIES 


1 . I Introduction 

Classical  multivariate  statistical  methods  deal  with  data  assuming 
values  in  Em,  Euclidean  space  of  m dimensions  consisting  of  m x 1 
vectors  with  real  components.  For  our  purposes,  we  consider  data 
restricted  to  Sm,  the  surface  of  the  unit  m-sphere  given  by 

sm  = [x.  : i £ Em,  x_'x_  = l}. 

The  case  m = l is  degenerate  in  as  much  as  = {—l , l } so  that  x e is 

a Bernoulli  random  variable.  For  m > 2,  Sm  is  the  domain  of  what  are 

commonly  called  "directional  data"  since  a vector  x = (x,,....x  )'  e S 

— i ’ ’ m m 

may  be  expressed  in  terms  of  the  "directions"  or  polar  angles  0l,02> 

* * * * ®m-l  by 


x 


i 


sin0^ 

sin0^ 

COS0j 


sin02* 

sin02* 


*sinem-l 

.sin0m_i  cos0m_.+l 


if  i = 1, 

if  2 < i < m- l, 

if  i = m, 


where  to  assure  a one-to-one  transformation  we  require  0 < 9^  < ir,  I < j 
< m-2  and  0 < 9m_1  < 2tt. 
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Most  instances  in  which  directional  data  arise  involve  observing 
the  motion  of  some  experimental  object.  Consequently,  much  attention 
has  historically  been  given  to  the  circular  case,  m = 2,  and  the 
spherical  case,  m = 3,  while  relatively  little  work,  has  been  devoted  to 
higher  dimensional  cases,  m > 3.  For  an  account  of  much  of  the  current 
statistical  methodology  surrounding  directional  data  on  the  circle  (2- 
sphere)  and  sphere  (3-sphere)  one  may  refer  to  Mardia  (1972). 

A statistical  methodology  dichotomized  according  to  directional 
data  on  the  circle  and  sphere  is  undesirable  for  two  reasons.  First, 
although  comparably  rare  at  present,  there  exist  applications  of 
directional  data  on  m-spheres  with  m > 3 as  indicated  in  section  1.2. 
Second,  there  is  a certain  redundancy  in  the  literature  in  which  methods 
now  treated  separately  for  the  circle  and  sphere  may  be  unified  under  a 
more  general  framework. 

For  example,  the  von  Mises  distribution  on  S£  and  the  Fisher 
distribution  on  are  treated  separately  in  the  literature.  However, 
both  are  special  cases  of  a general  distribution  on  Sm  introduced  by 
Watson  and  Williams  (1956)  and  now  called  the  Langevin  distribution  by 
Watson  (1983).  Moreover,  there  are  tests  in  the  literature  prescribed 
for  m = 2 and  tests  prescribed  for  m = 3 which  turn  out  to  be  exact 
analogues;  the  Rayleigh  test  is  one  example.  Elsewhere  there  are  test 
prescriptions  for  the  circle  for  which  there  are  no  immediately  apparent 
counterparts  for  the  sphere;  the  range  test  is  one  instance.  It  is 
revealing  that  tterdia's  text  treats  the  circular  and  spherical  cases 
separately.  While  readily  admitting  that  most  applications  occur  for 
m=2orm=3,  a desire  to  study  these  two  cases  simultaneously  and 
unify  certain  methods,  necessitates  the  consideration  of  procedures 
applicable  to  "general  m,"  m > 2. 


3 


1.2  Occurrence  of  Directional  Data 

Examples  of  directional  data  arise  from  a wide  variety  of 
sources.  In  astronomy,  directional  methods  are  used  to  study  the  motion 
of  planets  and  comets;  the  orbital  planes  of  our  sun's  planets  have  been 
the  subject  of  statistical  analysis  from  Daniel  Bernoulli  (1734)  to  G. 

S.  Watson  (1970).  Examples  in  biology  include  bird  migration,  swimming 
direction  of  daphnia  in  polarized  light  and  orientation  of  turtles  after 
treatment.  In  geography,  directional  data  are  often  recorded  in  terms 
of  longitude  and  latitude,  for  instance,  the  coordinates  of  earthquake 
epicenters.  Directional  data  are  used  in  geology  to  understand  remanent 
magnetism  (pole  wandering),  tectonic  forces  (such  as  continental  drift), 
lava  flow  patterns  and  orientation  of  cross  bedding  for  interpreting 
directions  of  depositing  currents  of  wind  or  water.  Of  interest  in 
medicine  is  vector  cardiographic  data,  monthly  onset  of  a disease  and 
monthly  death  rates.  Ifeteorological  studies  employ  directional  data  to 
describe  variables  such  as  wind  directions,  monthly  rainfall  and  monthly 
evaporation  from  a reservoir.  In  physics,  directional  data  have  been 
used  to  study  bubble  chamber  data,  atomic  weights  and  Brownian  motion. 

The  perception  of  direction  is  studied  in  physiology  and  psychology. 
Additionally,  circular  distributions  are  used  to  describe  such  diverse 
topics  as  certain  periodic  economic  data,  roulette  wheels,  rounding 
errors  in  numerical  calculations  and  the  generation  of  random  numbers 
(due  to  the  periodicity  of  remainders).  For  a discussion  of  the 
examples  above,  or  for  actual  data  sets,  one  may  refer  to  Batschelet 
(1965),  Mirdia  (1972)  and  Watson  (1983). 

On  occasion  we  find  directional  data  problems  which  require  more 
than  three  dimensions  for  analysis.  For  instance,  P.  J.  Laycock  (see 
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Mardia  (1975))  cites  multivariate  data  taken  on  fifteen  orthodontic 
planes  on  human  skulls.  A linear  model  indexed  by  a "direction"  in 
fifteen  dimensions  was  then  examined.  Further,  Stephens  (1982)  proposed 
an  analysis  of  variance  technique  employing  the  von  MLses  distribution 
to  analyze  a set  of  m proportions  of  a continuum  such  as  time  or 
volume.  Stephens'  procedure  requires  analysis  on  Sm  for  m > 2. 

1.3  The  Uniform  Distribution  on  the  Hypersphere 

Of  central  importance  in  developing  directional  data  methods  for 
use  on  general  hyperspheres  Sm,  m > 2,  is  testing  whether  or  not  the 
data  exhibit  some  "preferred  direction."  In  particular,  the  uniform 
distribution  on  Sm  provides  the  null  hypothesis  of  no  preferred 
direction.  In  section  1.5  we  will  more  fully  describe  the  role  of  the 
uniform  distribution  in  directional  data  and  in  section  1.6  we  will 
survey  methods  used  to  test  uniformity  on  Sm. 

In  order  to  define  the  uniform  distribution  and  other  commonly  used 
distributions  on  Sm  and,  indeed,  in  order  to  carry  out  integration  on 
Sm,  we  first  need  to  examine  the  topology  of  Sm. 

We  begin  by  defining  the  "cap," 

cm(i-’u)  " { 2L  : 2L  £ sm>  Jl'i.  > « }, 

to  be  an  open  set  on  Sffl.  In  other  words,  a set  GCS^  is  open  if  and 
only  if  for  each  x e G,  there  exist  = £(x)  and  u = u(x)  such  that 

— e ^rn( ?(x) , u(x_))CG.  Now  let  Am  be  the  collection  of  all  "caps,"  that 
is, 


u|  < 1 } 


Am  = { Cm(!’u)  : 1 £ Sm‘ 
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and  define  Bm  to  be  the  a-algebra  generated  by  Am. 

To  complete  the  formation  of  a probability  space  on  Sm,  a suitable 

measure  must  be  defined.  Although  any  non-negative  measure,  p,  on  the 

a-algebra  Bm  with  p(Sm)  = l will  suffice,  we  adopt  an  invariant  measure 

v = m/2  - l,  in  which  U^V^(A)  = U^V^(PA)  for  P an  orthogonal 

matrix  and  A e Bm.  Mjreover,  we  will  see  that  our  choice  of  measure 

conforms  to  the  usual  notion  of  arc  length  on  the  circle  S2  and 

surface  area  of  spherical  caps  on  the  sphere  Sg. 

Specifically,  let  be  a non-negative  measure  defined  on  the  a~ 

algebra  Bm  of  subsets  of  S with 
in  m 


U(v)(Cm(£,u))  = < 


C 1 (1-x2) 

J B(l/2, 


v-l/2 


dx 


vf  1/2) 


if  -»  < u < -1, 


if  -l  < u < 1 , 


if  l < u < ®, 


where  B(a,b)  is  the  beta  function  with  arguments  a and  b.  Note  that 
u^V^(sm)  = l and  that  the  invariance  property  of  requires  P(  x e 

Q 

m(_£,u))  to  be  independent  of  For  future  reference  we  may  also  write 


U(v)(Cm(£,u)) 


J 


u ( 1 2\v-1/2 

(1-XZ) 

_lB(l/2,  vfl/2) 


dX 


= l - (u) , say 
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To  see  the  connection  between  arc  length  in  the  circular  case, 
m = 2 (v  = 0),  and  the  measure  observe  that  the  arc  length  of 

C2(_§.,u)  is  2cos  lu.  Correspondingly,  the  probability  measure  gives 


dx 

B(l/2,  1/2) 
u 

which,  upon  setting  x = cos 9, 

-l 

Scos  u ^ 

de 

0 

2cos  ^u 
2 TT 

arc  length  of  C2(£,u) 
circumference  of 


For  the  spherical  case,  m = 3 (v  = 1/2),  the  surface  area  of  the 
spherical  cap  C^C^C.u)  is  2n(l  - u)  . The  probability  measure 
yields 


U(l/2>(C. 


- V — 

J B(l/2,  1 


dx 


) 


2 TT  ( l - u) 
4tt 


surface  area  of  C^(^,u) 
surface  area  of  S„ 
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With  the  probability  space  (Sm,  Bm,  U^v^)  now  at  our  disposal,  we 
may  define  a density  f(u)  of  Sm  with  respect  to  U^v^.  Let  _u  e Sm,  then 
f(u^)  necessarily  satisfies: 

(1)  f(u_)  > 0 almost  everywhere  (a.e.),  that  is,  for  every  _u  e Sffl 

except  on  a set  of  zero  measure; 

(2)  r f(u)  du(v)u)  = i. 

bm 

Now  let  u be  a vector  valued  random  variable  taking  values  on  Sm 
with  constant  density  f(_u)  = l . Then  we  shall  say  that  u_  is  uniformly 
distributed  on  Sm  with  respect  to  the  probability  measure  denoted 

_u_  - U(Sm).  It  should  be  noted  that  the  invariant  measure  im(*)  = 

otjjj  v^(  • ) is  frequently  used  in  the  literature.  Here,  am  = fx>x=i  Idx 

- 2nm^ /T(m/2)  is  the  surface  area  of  the  unit  m-sphere.  In  this  case, 
im(Sm)  = Ojh  so  that  the  uniform  density  with  respect  to  the  measure  im 

becomes  the  constant  l/c^.  For  the  most  part,  we  shall  adhere  to  the 

notationally  simpler  probability  measure  v = m/2  - l. 

We  now  present  a useful  theorem  regarding  the  uniform  distribution 

on  S . 
m 

Theorem  1.3.1:  Let  z_  have  an  nr-variate  normal  distribution  with  mean 

vector  0 and  covariance  matrix  I (denoted  z ~N(0,  I)),  and  put 

m — m — 7 nr  * r 

l /2  9 

r = II z II  = (z  1 z ) and  _u_  =_z/r.  Then  r has  a chi-square 
distribution  with  m degrees  of  freedom  (denoted  r^  ~ y^),  and  u ~ 
U(S  ) independently  of  r^. 

Proof : From  Saw  (1970,  1973)  we  find  the  Jacobian  of  the  transformation 

from  z_  to  (u_,r)  is  j{  z > u,r  } = r111-1 . Then  for  r > 0,  u e S , 

— — 1 1 — in’ 


8 


f (u , r ) = (2u)“m/2  r”"1  exp(-r2/2) 


(dr  dim) 


2 

and  for  r > 0,  u e S , 

7 — m 7 

(dr2  dl.) 
(dr2  dU(v)). 

Since  _u  e Sm,  we  have  u_  - U(Sm)  and  r2  - x£.  Further,  we  see  that  u_  and 
2 

r are  independent. 


f (u,r2) 


(r2)m/2  1 exp(-r2/2) 

2m~f*  r(m/2) 


m 


(r2)m/2  1 exp(-r2/2) 

2m^  r(m/2) 


l • 4 Other  Distributions  on  the  Hypersphere 

Although  the  uniform  distribution  will  occupy  much  of  our 
attention,  there  are  other  distributions  and  results  which  will  prove 
useful.  In  this  section  we  will  list  some  of  these  results  and  catalog 
several  distributions  for  later  reference. 

In  the  following,  a random  variable  Y will  have  the  beta 
distribution  if  the  probability  density  function  of  Y is  of  the  form, 


p(y)  = 


(y-a)P  1 (b-y)q_l 
B(p,q)  (b-a)P+q‘* 


(dy) , 


a < y < b, 


with  p > 0,  q > 0 and  this  will  be  denoted  Y - Bja  bj (p ,q)  . Also, 
unless  otherwise  indicated,  v = m/2  - l for  given  m. 
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Theorem  1.4.1:  Let  u - U(Sn)  and  _£  e Sffl.  Then  u'£  - B t_lf  (vH/2, 

v+1 /2) . 

Proof;  By  Theorem  1.3.1  we  may  write  u_'_£  - z_'jj/ 1 z_l  with  z_  - Nm(0_,  Ij,)  • 

Consider  a transformation  from  z_  to  Pz_ in  which  P is  an  orthogonal 

matrix  with  first  column,  (P)>p  equal  to  j,.  Then  without  loss  of 

2 2 1/2 

generality,  we  may  take  _£  = (l  ,0, . . . ,0) ' . Thus  u_'_£  - z^/(z|+  +zm) 

It  follows  that  Ou'J.)2  - z2/(zf+*“+z2)  -B[0>1](l/2,  \H-l/2).  Since 
z ^ / |2l  is  distributed  symmetrically  about  zero,  we  deduce  that  _u’_£  - 

B[-l  v+l/2>* 

Theorem  1.4.1  may  be  used  to  demonstrate  the  pertinence  of  our 
choice  of  measure  U(v).  With _u  and  as  in  the  preceding  theorem,  we 

have 

Cl  <*h)v'1/2  u-*)"-l/2 

P(  u'S  > t ) - \ 2v  dX 

J B(v+l/2,  v+l/2)  2 

C 1 (i-x2)^172  r(2vfi) 

= \ 2 — dx 

J r(vH/2)  T(  v+l/2)  2zv 

and  upon  using  the  Gauss  duplication  formula  for  gamma  functions, 

f 1 (i-x2)v-l/2  r(v+l) 

- \ dx 

J rU/2)  rCv+i/2) 

r 1 u-x2)v’l/2 

’ it dx 


B ( 1/2,  v+l/2) 
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= l - w^(t) 

where  W^v^(t)  is  as  defined  in  section  1.3. 

The  following  remark,  due  to  Saw  (1978),  will  be  of  frequent  use. 

Let  g( Xu)  satisfy 

(1)  g( Xu)  > 0 a.e.  with  respect  to  the  measure  W^v\u) 

where  -l  < u < l; 

(2)  Jl  g( Xu)  dW(v)(u)  = l. 

Then  following  Saw  (1978),  we  show  that  g(Xu_'_£)  is  a density  with 
respect  to  U^v\u_).  Using  Theorems  1.3.1  and  1.4.1  we  have 

Is  gUi'i.)  dU(v)(u_)  = Eu(g(Xu_’0),  u - U(SB) 

m — 

= Ez(g(  Xz ’_§/  tlz_i ) ) , z_  - Nm( 0_,  Iffl) 

= Et(g(Xt)),  t - B[.lflJ(vfl/2,vH/2) 


-s 


v-l/2 


-l 


g(Xt)  (l-t2) 

B (1/2,  vfl/2) 


— dt 


= l. 


Hence  g(Xu_'_£)  is  also  a density  on  Sm  with  respect  to 
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Theorem  1 .4.2 : Let  _u  - U(Sm)  and  let  Q be  an  m x m symmetric  idempotent 
matrix  of  rank  k.  Then  u_'Qu  ~ Bj-q  (k/2 , (m-k) /2) . 

Proof : Choose  the  orthogonal  matrix  P such  that 


PQP' 


0 

0 


and  by  Theorem  1.3.1  write  u_  ~ _z_  / llzjl  where  _z_  - Nm(0_,Im).  Then  u_'Qu_  - 
(Pz)'Q(Pz)  / llzjl2  = z_'P'QP_z  / (z_'_z)  = (zJ+***+z£)  / (z2+**’+z2)  - 
B [0 , l ] (k/2  > (m-k) /2) . 

Theorem  1.4.3:  Let  _u_  - U(Sm)  and  _£  e Sm.  Then  (m-1)^2  u_'_£  / 

(l“(u_'jj)2)  ^2  has  a Student's  t-distribution  with  m-l  degrees  of 

freedom  (denoted  t„_,). 

m-  i 

Proof : Again  we  use  Theorem  1.3.1  and  as  in  the  proof  of  Theorem  1.4.1 

take  = (l,0,...,0)'  without  loss  of  generality.  Hence  we  have, 

(nr-l)l/2  uV^/Cl-CuVg)2)1/2  = z ^ / [ (z|+*  * *+z2  ) / (m- 1 ) ] 1 ^2  ~ ^n-l* 

In  future  work  we  will  also  have  occasion  to  appeal  to  a central 

limit  theorem  for  vector  valued  random  variables  which  we  record  here 
(see  Rao  (1973)).  Let  ^,...,2^  be  independent  identically  distributed 
(i.i.d.)  vectors  in  Em  with  Ej^  = jj  and  CovG^)  = E,  l < i < n.  If  _s  = 
x_l+...+xn,  then  (s_-np)  n_1/^2  converges  in  distribution  (^)  to  Nm(_0 , E) , 
that  is,  U-np)  n-1^2  ^ 

In  subsequent  work  we  shall  restrict  our  attention  to  those 
densities  on  Sm  which  are  of  the  form  f(Xu'_£)  and  which  satisfy 
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(1)  f(Au^_£)  is  constant  on  u_'j^  = t,  that  is,  f(Au^j^) 

. . . is  rotationally  symmetric  about 

v l *4 • l J 

(2)  f(Aw)  is  a non-decreasing  function  of  w over  [-1,1]. 

This  restriction  seems  appropriate  since  many  commonly  occurring 
densities  arising  in  practice  manifest  the  conditions  (1.4.1).  In  fact, 
the  conditions  (1.4.1)  guarantee  that  f(Au_'_£)  is  unimodal  with  mode 
located  at 

We  now  introduce  several  important  distributions  in  directional 

data. 


1.4.1  Langevin  Distribution 

The  Langevin  density  on  Sm  is  given  by 
U*4*2)  r(v)(Au'£)  = exp(  Au_'_£)  / Ay(  A)  (dU(v)(u)) 

where 


AvCA) 


\ 


1 a-,2,-1/2 

exp(Ax)  dx 

_x  B(l/2,  vfl/2) 


- ( A/2) v r(v+D  IV(A) 


and  IV(A)  is  the  modified  Bessel  function  of  order  v given  in  series 


form  as 


u - (A/2)2k 

I (A)  = (A/2) v l . 

k=0  k!  r(v+k+l) 
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From  this  we  see  that  Ay(A)  may  also  be  expressed  in  terms  of 

2 

hypergeometric  functions  as  Ay(A)  = qFjC-jvHjA  /4).  The  quantities^ 
and  l /A  are  somewhat  analogous  to  the  ordinary  mean  and  variance  on  the 
real  line. 

The  cases  m = 2 and  m = 3 are  of  special  interest.  On  the  circle, 
m = 2,  the  Langevin  distribution  is  commonly  called  the  von  MLses 
distribution  owing  to  the  fact  that  von  Mises  (1918)  introduced 
r^^(Au_'Jj)  = exp(Au^_g)  / Iq(A)  as  a density  on  the  circle  in  order  to 
study  deviations  from  integral  values  of  measured  atomic  weights. 

Since  significant  advances  for  the  case  m = 3 were  made  by  Fisher 
(1953)  while  working  on  a paleomagnetism  problem,  the  Langevin 
distribution  on  the  sphere  is  known  as  the  Fisher  distribution  with 
density  r(  1 /2) ( Au/£)  = (A/sinhA)  exp(Au_'_£).  The  density  r(v)(Au/j;)  was 
generalized  to  Sm  by  Watson  and  Williams  (1956).  But  the  first  to 
examine  the  density  r^ 1 ^ ^ ( Au_'_£)  seems  to  have  been  Langevin  (1905)  in 
regards  to  a statistical  mechanics  study  of  magnetism.  Hence,  Watson 
(1983)  has  proposed  the  name  "Langevin"  distribution  to  settle  any 
question  of  authorship  priority. 

If  u_  has  density  r( v) ( Au_'_§)  we  write  u_  - R(v)(A,_£).  A few  notable 
properties  of  the  Langevin  distribution  are  as  follows: 

(1)  for  s = xl^~***+xn  sum  n unit  vectors  on  Sm,  then 

A* 

S.  = sfist  is  the  maximum  likelihood  estimator  of  if  and  only 
if  - i.i.d.  R(v)(A,j;),  1 < i < n; 

(2)  as  A 0,  R^V^(A,^)  tends  to  the  uniform  distribution  on  S ; 

m’ 
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(3)  the  von  Mises  distribution  on  the  circle  maximizes  entropy 
under  the  constraint  that  the  mean  direction  and  circular 
variance  are  fixed. 

It  should  be  noted  that  the  Langevin  distribution  shares  many 
properties  analogous  to  the  normal  distribution  on  the  real  line.  For 
example,  properties  (1)  and  (3)  above  have  analogues  in  the  linear 
normal  case. 

1.4.2  Saw's  Distribution 

A density  on  Sm  introduced  by  Saw  (1978)  is  defined  by 


(1.4.3)  j^V^(Xu'£)  = exp(-X2)  1 


(2Xu'£)k  r(v+k/2+l) 


k=0 


k!  r(v+l) 


If  _u  has  density  j^(Xu_'j,)  we  write  _u  - J^(X,_£). 


Many  interesting  properties  of  the  distribution  J^(X,_£)  are  given 
by  Saw  (1978).  Among  these  properties  are  the  following: 


(2)  as  X + 0,  J^v\x,S)  tends  to  the  uniform  distribution  on  S : 

— m’ 


(3)  if  _x  - Nm(JL>  ^ and  _U  t 0_,  then  u_  = x/llxjl  - J^(X,_£)  where 
= jj/ HjjJI  and  X = Il2_^^_ull. 
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1.4.3  Saw's  Hyperspherical  t-Distribution 

Saw  (1984)  proposed  the  density  on  Sffl, 

(1.4.4)  e^v)(Au'l)  =KS(X)  / (l-2Au'i+A2)X+s  (dU(v)) 

for  s = 1,2,...,  where  in  terms  of  hypergeometric  functions, 

k“1(\)  = 2f  i (s » v+s ; v*-1 ; ) 

= (l-X2)-s  21" i Cs , l— s ; \rt-l ; X^/ Cl — ) ) . 

If  u_has  the  density  e^(Au_'_£)  we  write  u_  - E^(A,_£). 

Some  properties  discussed  by  Saw  (1984)  include 

(1)  if  uj_£  - E^v)(A,£)  and  j,  - e(v^(p,jt),  then  u has  the  Cauchy- 
like  property  that  u.  - E|V^(Ap,jt); 

(2)  as  X + 0,  EgV^ ( A,jj)  tends  to  the  uniform  distribution  on  Sm; 

(3)  as  s * ®,  E^(A/s,_£)  tends  to  the  Langevin  distribution 

R(v)(X,l). 

We  note  that  e[v^(X,_£)  is  the  generalization  of  the  Cauchy 
distribution  to  Sffl;  here  e^v^(Xu_'_£)  = (1~A2)/ (l-2Xu^jj+A2) . In  fact 
for  s = 1 and  v = 0,  e[°\a,_£)  is  the  circular  wrapped  Cauchy 
distribution  (or  essentially  the  Poisson  Kernel  from  complex  analysis). 
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Now  consider  the  family  of  densities 


T = { egV^(Xu_'jj)  : X > 0,  e Sm,  s = 0,1,2,...  } 


which  contains  the  Cauchy  distribution  on  one  end  of  the  spectrum  for  s 
and  the  Langevin  distribution  (analogous  to  the  linear  normal)  on  the 
other  end.  Thus  we  see  that  the  class  T shares  a property  similar  to 
the  Student's  t-distribution  in  that  T contains  distributions  analogous 
to  the  linear  Cauchy  and  linear  normal  distributions.  It  is  in  this 
sense  that  we  call  a distribution  from  the  class  T a "hyperspherical  t- 
distribution. " 

The  following  three  distributions  are  given  only  for  the  circular 
case.  Due  to  the  need  for  Gegenbauer  polymonials  to  express  these 
densities  on  Sm,  we  postpone  the  generalizations  to  Sm  until  Chapter 
2.  Rather  than  describe  these  circular  distributions  in  Cartesian 
coordinates,  x_  = (x^,X2)'  e S2,  we  use  the  conventional  polar  angle 
notation  for  circular  distributions,  9 e ( 0 , 2 n ] , where 


r 


x^  >0,  *2  > 0, 


ir/2, 


xL  = 0,  x2  = l, 


3ir/2, 


xt  =0,  x2  = -1, 


and  -ir/2  < tan-1  (x2/xp  < ir/2. 
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1.4.4  Circular  Cardioid  Distribution 

The  circular  cardioid  distribution,  introduced  by  Jeffreys  (1948), 
has  density 

(1.4.6)  f c( 0)  = l + 2Xcos(6-e0)  (dU(0)) 

where  0 < 9 < 2tt,  0 < 9q  < 2tt  and  | X | < 1/2.  Observe  that  for  \ = 0 
the  cardioid  distribution  reduces  to  the  uniform  distribution  on  $2* 

1.4.5  Circular  Point  Distribution 

The  density  on  S2  which  concentrates  all  the  mass  at  the  point 
9 = 9q  is 

00 

(1.4.5)  f (9)  =1+2  l cos  p(9-9„)  (dU(0)) 

P p-l 

where  0 < 9 < tv. 

1.4.5  Circular  Wrapped  Normal  Distribution 

The  circular  wrapped  normal  distribution  has  density 

,/o  “ “I 

(1.4.7)  f (9)  = (2ir/o2)  £ exp{ ( 9-9n+2kir)2} 

W k=-»  2a2  U 

with  a > 0,  0 < 9q  < 2ir  and  0 < 9 < 2 tt . We  may  also  express 


Cdu<0)) 
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f (e)  = t + 2 l Pk  cos  k(e-en)  (du(0)) 

W k=l 

where  p = exp(-02/2)  and  0 < 9q  < 2 tt . 

Among  properties  for  the  circular  wrapped  normal  distribution  are 

(1)  as  a * 00  (p  -►  0),  f wC 0 ) tends  to  the  uniform  density  on  S2 ; 

(2)  as  a 0 ( p * l)  , fw(9)  tends  to  the  circular  point 
distribution  on  9 = 9q ; 

(3)  if  9j,...,9  are  independent  wrapped  normal  random  variables, 

each  with  mode  0 and  respective  parameters  then 

9^+.  • *+9n  (mod  2 tt ) has  a wrapped  normal  distribution  with 

2 2 2 

parameter  a = a^+* • and  mode  0 ; 

(4)  the  circular  wrapped  normal  distribution  provides  the 
distribution  for  Brownian  motion  on  the  circle; 

(5)  if  9j  - i.i.d.,  j =l,...,n,  with  common  distribution  function 

F(9),  E(9)  = 0 and  — ir  < 9.  < tt  for  each  i,  then  s = 

J J ’ n 

1/2 

n (8j+***+9n)  (mod  2 tt ) tends  to  the  wrapped  normal 
distribution  with  a ^ = E(92). 
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1.5  The  Role  of  the  Uniform  Distribution 

As  already  remarked  in  section  1.3,  the  uniform  distribution 
provides  the  null  hypothesis  of  no  preferred  direction  for  data  on  Sm. 

It  is  in  this  setting  that  we  are  interested  in  testing  for  uniformity 
on  S . The  most  prevalent  prescriptions  for  testing  uniformity  on  Sm 
are  discussed  in  section  1.6. 

But  beyond  being  of  use  in  providing  the  null  hypothesis  of  no 

preferred  direction,  the  uniform  distribution  arises  in  a variety  of 

ways  in  directional  data.  Indeed  the  uniform  distribution  shares  some 

properties  analogous  to  those  of  the  normal  distribution  on  the  real 

line,  albeit  most  of  the  important  analogous  properties  are  possessed  by 

the  Langevin  and  wrapped  normal  distributions.  In  this  section  we  shall 

examine  some  characterizations  of  the  uniform  distribution. 

One  property  of  the  normal  distribution  on  the  real  line  due  to 

Geary  (1936)  is  that  for  n > 2 i.i.d.  observations,  the  mean  and 

variance  are  independent  if  and  only  if  the  observations  are  normally 

distributed.  A similar  result  holds  for  observations  on  S„. 

m 

l* t x.j , . . . ,2^ , (n  > 2),  be  independent  observations  on  Sm  with 
common  density  f(x_)  which  is  continuous  almost  everywhere  (a.e.)  with 
respect  to  Lebesgue  measure  on  Sm.  Wit h _s_  = 2Li+*  * *+2En  > define  the 
direction,  _Xq,  of  the  sample  mean  and  the  resultant  length,  rn,  by 

(1.5.1)  r = II  s II  and  xn  = s / r . 

The  quantities  _xq  and  n-rn  are  somewhat  analogous  to  the  sample  mean  and 
variance  on  the  real  line. 

With  Lebesgue  measure  on  Sm,  Bm  as  defined  in  section  1.3  and 
representing  the  surface  area  of  Sm,  let  b = ess  inffa^  f(_x)  : _x  e Sm}  = 
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sup(inf{am  f(x_)  : x_  e Sm})  where  the  supremum  is  taken  over  the  class  of 

all  real  valued  B measurable  functions  on  S which  are  equal  to  f(x) 
m m 

a.e.  and  the  infimum  is  taken  over  all  _x  e Sm.  Then  Kent,  Mardia  and 
Rao  (1979)  have  shown  that  among  all  distributions  with  a density  and 
for  b > 0,  Xq  and  rQ  are  independent  if  and  only  if  f(x)  is  the  uniform 
density  on  Sm.  The  rather  artificial  condition  that  b > 0 is  also  shown 
by  Kent  et  al.  (1979)  to  be  unnecessary  for  ra  = 2. 

The  uniform  distribution  analogously  shares  another  property  with 
the  normal  distribution  on  the  real  line  for  it  appears  as  the  limiting 
distribution  in  the  central  limit  theorem  on  the  circle. 

More  precisely,  let  the  angles  t>e  i*i*d.  random  variables 

with  common  circular  distribution  function  F(0).  Then,  provided  F(9)  is 
not  the  lattice  distribution  with  mass  distributed  on  the  vertices  of  a 
regular  polygon  inscribed  in  the  unit  circle  with  no  lattice  point 
coinciding  to  zero,  the  sum  0^+»**+0n  (mod  2 tt ) converges  in  distribution 
to  the  uniform  distribution  on  the  circle.  Also  under  the  constraint 
that  the  random  variables  possess  axial  symmetry,  that  is,  f(_x)  = f(-_x), 
_x_  e S^,  the  central  limit  theorem  on  the  sphere  gives  the  sum  of  random 
variables  converging  to  the  uniform  distribution  (see  Breitenberger 
(1963)). 


1 . 6 Testing  Uniformity 

A fundamental  inference  problem  for  data  on  Sm  is  testing  whether 
or  not  the  data  exhibit  a preferred  direction,  that  is,  determining 
whether  or  not  the  data  originate  from  a uniform  distribution.  Besides 
having  practical  importance,  testing  uniformity  is  of  central 


theoretical  interest. 
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To  accommodate  applications  in  which  higher  dimensional  spheres  are 
necessary  and  more  essentially  to  provide  a unified  methodology,  we 
examine  tests  of  uniformity  which  encompass  two,  three  or  more 
dimensions.  The  bulk  of  proposed  tests  of  uniformity  have  been 
presented  in  two  or  three  dimensions.  Indeed,  Mardia  (1975)  notes,  "As 
in  multivariate  analysis,  it  is  not  surprising  that  there  is  a dearth  of 
nonparametric  tests  of  practical  value  for  p > 3”  (p.  369).  In  this 
section  we  present  some  conventional  tests  of  uniformity  in  two  or  three 
dimensions  which  have  extensions  to  higher  dimensions,  m > 3. 

Historically,  the  earliest  test  of  uniformity  was  based  on  the 
vector  resultant.  The  distribution  of  the  resultant  was  studied  by  Lord 
Rayleigh  (1880)  in  regards  to  a problem  in  sound  concerning  the 
composition  of  n iso-periodic  vibrations  of  unit  amplitude  and  random 
phases . 

Apparently  Rayleigh's  solution  was  not  widely  known  among 
statisticians  when  Karl  Pearson  (1905)  posed  the  problem  of  the 
distribution  of  the  resultant  length  of  a man  taking  n equal  length 
steps  from  the  origin  where  the  direction  of  each  step  is  chosen  at 
random.  Upon  learning  of  Rayleigh's  solution  for  large  n,  Pearson 
(1905)  colorfully  remarked,  "The  lesson  of  Lord  Rayleigh's  solution  is 
that  in  open  country  the  most  probable  place  to  find  a drunken  man  who 
is  at  all  capable  of  keeping  on  his  feet  is  somewhere  near  his  starting 
point!”  (p.  342). 

Rayleigh's  test  of  uniformity  in  m dimensions  may  be  based  on  an 
approximation  to  the  distribution  of  the  resultant  when  the  sample  size, 
n,  is  large  (see  Watson  (1983)  for  instance). 

Suppose  x^, ... ,x^  are  i.i.d.  observations  taken  from  the  uniform 
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distribution  on  Sm»  Let  _p_  - Ex  and  Z E (x~y ) (x~p ) • Since  P_x_  - U(Sm) 

for  any  orthogonal  matrix  P,  we  have  jji  = P_y  so  that  _p_  = _0_.  In  addition, 

since  Z = PZP'  and  trace(Z)  = trace(Exx'  ) = E trace(x.'x)  = 1 it  follows 

that  Z = Im/m.  Now  withj^=  x^+« » «+xn  and  using  the  central  limit 

theorem  for  unit  vector  random  variables  listed  in  section  1.4,  we  find 

_s/nl/2  i Nm(0>Vm)*  Thus  (m/n)l/2  s_  * Nm(  0_,  Im)  and  putting  rn  = 

(s's)1/2,  we  may  write  (m/n)  r2  as  the  test  statistic  and  reject  for 

2 

large  values  according  to  the  percentage  points  of  the  xm  distribution. 

Among  tests  for  uniformity  on  the  circle  which  have  received  much 
attention  are  those  due  to  B.  Ajne  (1968).  Suppose  n points  are  sampled 
independently  from  the  circle  S2  with  Lebesgue  measure.  Arbitrarily 
select  an  origin  on  S2,  say  Q,  and  select  a positive  direction,  say 
counterclockwise.  An  arc  on  S2  of  length  0 may  be  expressed  as  those 
points  on  S2  whose  distance  from  Q belong  to  the  interval  (x,x+9] , 0 < x 
< 2tt . Define  N(9,x)  to  be  the  number  of  observations  belonging  to  the 
arc  corresponding  to  the  interval  (x,x+9] . For  example,  in  Figure  l.l 
with  n = 7 sample  points,  we  see  that  N(0,x)  = 4. 


Figure  l.l  Sample  Points  in  an  Arc 


of  Length  9 on  the  Circle 
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Now  set 


(1.6.1) 


A(  0)  = 


2im 


S 


2tt 


( N ( 9 , x ) 


n9  s2 
" ~2-n~  > 


dx 


and 


(t .6.2) 


B(0)  = max  N(9,x). 
0<x<2ir 


Ajne  (1968)  studied  the  statistics  I = 2irn  A(ir)  and  N = B(ir)  while 
Rothman  (1971)  presented  slightly  more  general  results  for  A(9)  and  B(9) 
with  0 < 9 < 2tt.  Further  generalizations  have  been  given  by  Be  ran 
(1968,  1969),  Gine  (1975),  Prentice  (1978),  Jupp  and  Spurr  (1983)  and 
Saw  ( l 984) . 

Ajne  was  concerned  with  testing  the  null  hypothesis  of  uniformity 

1 

Hq  : f(0)  0 < 0 < 2tt 

2 TT 


against  the  alternative 


Hj : f ( 0) 


0 < 9 < it 
it  < 0 < 2ir 


where  p + q = l and  0 < p < 1. 

Ajne  found  that  rejecting  Hq  for  large  values  of  N = B(ir)  provides 
the  most  powerful  invariant  (under  rotations)  test  against  "distant" 
alternatives  in  H^  with  p near  1.  Using  a random  walk  argument,  Ajne 
found  the  probability  distribution  of  N to  be 
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P(  N > k)  = 2_(n  (2k  - n) 


j (2k  - a)  + k 


a 


where  k = [n/2]  + and  [ • ] is  Che  greatest  integer  function. 

These  probabilities  were  tabulated  by  Ajne  (1968)  and  also  by  Klotz 
(1959)  as  probabilities  associated  with  a bivariate  sign  test  offered  by 
Hodges  (1955).  The  connection  between  Aj ne' s test  and  Hodge's  test  was 
first  noted  by  Bhattacharyya  and  Johnson  (1969). 

The  asymptotic  distribution  of  N = (2N  — n)  / n as  the  sample 
size  n tends  to  infinity  was  given  by  Ajne  (1968)  as 


where  L(z)  is  the  limiting  distribution  of  the  Kolmogorov-Smirnov 
statistic  D(n),  (see  Feller  (1948)).  The  power  function  for  N and 
asymptotic  power  function  for  N were  computed  and  recorded  by  Ajne 
(1968). 

Of  the  two  statistics  N and  I introduced  by  Ajne,  it  is  the  variant 
A(n)  = 1/(2™)  which  has  been  the  subject  of  much  generalization.  For 
alternatives  in  H^  near  uniformity,  that  is,  as  p * 1/2  in  Hj,  Ajne 
(1968)  established  that  rejecting  Hq  for  large  values  of  I gives  the 
locally  most  powerful  invariant  (under  rotations)  test. 

Using  a Fourier  expansion  technique  proposed  by  Watson  (1967),  the 
limiting  distribution  of  2I/(ira)  as  n * <*>  is 


lim  P(N*  > c)  = L(tt/(2c)) 


“ L 
7T  2 m=l 


(2m  - l)2 


— 1/2  2 
in  which  p = 1/2  + xn  ' , x > 0 and  the  random  variables  ara 
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independent  chi-square  variates  with  two  degrees  of  freedom  and  non- 

2 2 2 

centrality  parameter  X = 16x  / [(2m-  l)  ir  ].  This  limiting 

2 2 

distribution  has  mean  l + 4x  and  variance  2 + I6x  /3.  Under 
uniformity;  that  is,  when  x = 0,  the  limiting  distribution  assumes  the 
form 


P( 


21 

nir 


> c ) 


(-!)■*  exp [-(2j+l ) 2 tt2c/8] 

23  + 1 


c > 0. 


Following  Prentice  (1978),  the  m-dimensional  analogue  of  Ajne's 
statistic  I/(2irn)  for  a random  sample  _xp  •••  fr°m  Sm  (with  Lebesgue 
measure)  is 

(U6*3)  Am,n  = (nam)_l  / , j (N(l)-n/2)2  di^) 

where  am  = 2irm/^2  / T(m/2)  is  the  surface  area  of  the  unit  m-sphere  and 
N(_£)  is  the  number  of  data  points  x^  such  that  Jj/iLi  > 0.  Alternately  we 
may  express 


A 

m,n 


(rm)  V V cos  ^(x!x.)  . 

i<j  -‘-J 


The  asymptotic  null  distribution  of  A is 

m,n 


00 


l 

q=l 


2 2 

a2q-l  Xd(2v+l,2q-l) 


where 


a 


2q— 1 


(_l)q  1 22v  r(vfl)  T(q+v)  (2q-2) ! 


ir  (q-l)  ! (2q  + 2v  - 1)  ! 
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and  the  X2d(2v+1  2q-l)  random  variables  are  independent  with  degrees  of 
freedom 


d(2vH , 2,-D  - (2q+^2)  + 


Beran  (1968,  1969)  considered  generalizations  of  Ajne's  statistic 
A(tt)  by  proposing  a class  of  test  statistics  on  the  circle  S2 

1 f 1 n 2 

(1.6.4)  T * \ [ l f(*4-x  ) - n ] dx 

n n J 0 j=l  J 

where  f(x)  is  a non- uniform  probability  density  on  S2  (with  respect  to 
the  measure  U^).  Beran  showed  that  in  special  cases,  Tn  reduces  to 
Rayleigh's  test  and  Ajne's  test  among  others. 

Both  Saw  (1984)  and  Prentice  (1978)  have  studied  generalizations  of 
the  statistics  Tn  to  m-dimensions.  Saw  (1984)  gives  this  generalization 
as 


(1.6.5)  G(v)  = [ II  4>(xU)  - 1 ]2  dU(v)U) 

n J S j=l  J 

where  _Xj  , ...  ,x^  is  a random  sample  with  common  density  <}>(x_)  over  Sm 
(with  respect  to  the  measure  U^v^).  For  certain  choices  of  ^(jc),  the 
statistics  G^v^  can  be  shown  to  include  Rayleigh's  test  and  Ajne's  test 
among  others.  Using  Gegenbauer  polynomials  (see  Chapter  2),  it  is 
possible  to  express  G^v^  in  summation  form.  According  to  Saw  (1984), 


,(v) 


n 

I 

i=l 


n 

I 

j = l 


I 

p = l 


( >2 


» ( V ) f , V 

A (u.u.) 
P -i“J 
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v + p [ p/2 J (-l)S  r(v+p-s)  (2x)P_2s 

£ , p * o, 

r(\H-l)  s=0  s!  (p-2s)  ! 

I, 

A(v)(l)  -Av)  = j +(_u^)  Ap(v)(u'l)  dU(v)(i_). 

P P Sm 

Prentice  (1978)  also  gives  the  asymptotic  distribution  of  under  the 

null  hypothesis  of  uniformity. 

Another  class  of  tests  which,  for  the  cases  m = 2 and  m = 3, 
includes  the  tests  of  Rayleigh,  Ajne,  Beran  and  others  has  been  given  by 
Gine  (1975).  Gine's  class  of  tests  is  based  on  Sobolev  norms  and  has 
since  been  generalized  further  by  Jupp  and  Spurr  (1983)  who  considered 
tests  for  symmetry  also  based  on  Sobolev  norms. 


in  which 


A(U)U)  - 
P 


A<v)(x)  = 


CHAPTER  2 

GEGENBAUER  POLYNOMIALS  AND  DENSITY  REPRESENTATION 
ON  THE  HYPERSPHERE 

2.1  Introduction 

The  representation  of  densities  on  Sm  will  proceed  by  use  of 
spherical  harmonic  functions.  Ihus  we  present  a brief  description  of 
some  important  notions  from  spherical  harmonics. 

Let  Gn(x) , x.  £ Em  be  a homogeneous  polynomial  of  degree  n in  m 
dimensions.  That  is, 

Gn(rx)  = tn  Gn(x), 

x e Em  and  t a non-zero  scalar.  If  Gn(_x)  satisfies 

\ =„<*>  - 0 

where 

7»  ‘ M (’k)  ' (-^r)  + " + ("%“) 

is  the  Laplacian  operator,  then  Gn(_x^)  is  called  a harmonic  function. 

With  u e S , define 
— m 

Sn(m,u)  = Gn(u). 

Muller  (1966)  calls  Sn(m,_u),  ju  e Sm,  a spherical  harmonic  function  of 
order  n in  m dimensions.  The  spherical  harmonic  functions  are 
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orthogonal  with  respect  to  the  measure  im  in  that,  for  e Sm  and  p ¥ q, 


/ S (m,u)  S (m,u)  di  (u)  = 0. 

S P — 9 — m — 

m 

In  addition,  there  exist  N(m,n)  = [(2n+m-2)  r(n+m-2) ] / [r(n+l)  r(m-l)] 
linearly  independent  spherical  harmonics  S . (m,u),  j = l , . . . ,N(m,n) , 

n » j 

forming  a basis  for  the  spherical  harmonic  functions  of  degree  n in  m 

dimensions.  We  may  suppose  that  S . (m,u),  j = l,...,N(m,n)  are 

^ > J 

orthonormal;  that  is, 

/s  SpJ(m,u_)  Sq  ,j  ^m»— ^ dim^  = 6pq 
in 

where  the  Kronecker  delta,  6p^,  is  defined  to  be  zero  for  p ^ q and 
unity  for  p = q.  Thus  for  any  _u  e Sm  we  may  write 

N(m,n) 

S (m,u)  = V c . S ,(m,u) 
n " n,J  n,  J “ 


for  some  scalars  (c  .}  determined  by 
1 n,jJ  J 

cn,  j = /s  Sn,j(m‘ii.)  • 

m 

However,  the  basis  elements  Sn  j(m,ja^)  are  complicated  to  represent  for 
general  m-dimensional  spheres.  To  achieve  simpler  representations  we 
restrict  attention  to  the  subset  of  homogeneous  harmonic  polynomials  of 
degree  n in  m dimensions,  Pn(m,x_),  x_  e Em,  with  the  property  that 


Pn(m,Ax)  = Pn(m,jc) 
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for  all  orthogonal  transformations  A which  leave  one  vector,  say  _z_, 
fixed.  The  polynomials  Pn(m,_z)  are  unique  up  to  a multiplicative 
constant.  The  constant  may  be  fixed  by  requiring  Pn(m,jO  = 1 where  the 
vector  z_  is  invariant  under  the  orthogonal  transformations  described 
above.  Muller  (1966)  terms  {PnCm,x_) } Legendre  functions  of  degree  n in 
m dimensions  and  shows  that  every  spherical  harmonic  function  Sn(m,^j_), 
u_  e Sm,  may  be  written 

N(m,n) 

S (m,u)  = T a.  P (m,  n!u) 
n - J n -j- 

where  rh,...,n  are  unit  vectors  with 

~TJ(m,n) 

det{Pn(m,_n;[_nj)  : i,j  = 1 , . . . ,N(m,n)  } > 0. 

We  shall  be  concerned  with  the  expansion  of  functions  in  terms  of 
spherical  harmonic  functions.  It  is  shown  by  Muller  (1966)  that  every 
continuous  function  F(ji)  , _u_  e Sffl,  can  be  approximated  uniformly  by 
spherical  harmonic  functions  in  that, 

00 

F(u)  = lim  y xn  S (m,u) 

T+l-0  n=0  n 


where 


sn(m,u_)  = c^1  N(m,n)  / Pn(m,_u'_£)  F(_g_)  difflU) 

^m 

and  is  the  surface  area-  of  Sm.  Further,  Miller  (1966)  shows  that  the 
system  of  spherical  harmonics  is  complete  for  continuous  functions  on  S 


in  the  sense  that  if 
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/ F(u)  S (m,u)  di  (u)  = 0 

for  all  spherical  harmonics  Sn(m,_u),  then  F(_u_)  is  identically  zero. 

Fortunately,  it  suffices  for  our  present  purpose,  to  restrict 
attention  to  the  subset  (Sn(m,_£)}  having  the  property  that  Sn(m,A_g)  = 
Sn(m,_£)  for  all  orthogonal  matrices  A such  that  Ajj^  = It  will  also  be 
important  that  e Sm,  can  be  expressed  in  terms  of  the 

Gegenbauer  polynomials  defined  in  the  next  section.  Through  the  use  of 
Gegenbauer  polynomials  we  will  achieve  representations  for  densities  of 
the  form  fdu.'jj.)  which  are  restricted  to  be  rotationally  symmetric  about 
e Sm;  that  is,  f(Xu/_£)  is  invariant  under  orthogonal  transformations 
which  leave  one  point,  jj_,  on  the  m-sphere  fixed.  In  addition  we  require 
to  be  nondecreasing  in  £'^.  But  before  providing  such  a density 
decomposition  in  section  2.3,  we  examine  many  of  the  standard  properties 
of  the  Gegenbauer  polynomials  as  collected  in  the  next  section. 

2.2  Gegenbauer  Polynomials 

We  define  {C^(x)J,  n = 0,1,2,...  to  be  a system  of  polynomials 
orthogonal  on  the  interval  [-1,1]  with  respect  to  the  weight  function 
(l-x2)v_1/2,  v > -1/2;  that  is,  for  m ^ n. 


The  weight  function  (l-*2)v-l/2  d.tetlIlMS  the  poly„omials  (cu(x) j apart 

from  a constant  factor.  Commonly  called  Gegenbauer  or  ultraspherical 
polynomials,  the  polynomials  {C^(x) } are  part  of  the  well  known 
classical  orthogonal  polynomials  along  with  Hermite,  Laguerre,  Legendre 
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and  Jacobi  polynomials.  The  notation  "C^(x)"  is  due  to  Leopold 
Gegenbauer  who  studied  these  polynomials  in  a series  of  papers  from  1877 
to  1893. 

Being  unique  only  up  to  an  arbitrary  multiplicative  constant, 
orthogonal  polynomials  are  usually  standardized  in  some  way.  The 
conventional  standardization  is  to  put 

/ n + 2 v - 1 \ 

cJJ(l)  - I j > v * 0 

C°(l)  = , n * 0 

n a 

cgd)  - 1. 

However,  we  choose  a standardization  given  by  Saw  (1984)  which 
lends  itself  to  simplifying  certain  properties  of  Gegenbauer 
polynomials.  For  example.  Saw's  choice  of  standardization  and  weight 
function  allows  the  integrated  square  of  a Gegenbauer  polynomial  to 
assume  the  particularly  simple  form  given  in  (2.2.67). 

Saw  (1984)  defines  the  polynomials  (ApV\x)}  by 

r v + p L p/2  j (-i)s  r(v+p-s)  (2x)p  2s 

A V (x)  l , p * 0, 

P r(v+l)  s=0  s!  (p-2s)! 

A^v)(l)  = 1. 

Equivalently  we  may  define  the  polynomials  |a^v^(x)}  by  the  generating 
relation  given  by  Saw  (1984)  as 


I - X2 

(l  - 2Ax  + X*)*T 
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OO 

= y A(v)(x)  XP  , 111  < l. 

P=0  p 11 

The  polynomials  (a£v\x)}  are  related  to  {Cp(x)}  by 

( p+v 

A^V;(x)  CV(x)  , v 4 0, 

P v P 

= P Cp(x)  , P 4 o, 

= 1 , v = p = 0. 

Beyond  simplifying  certain  results,  the  standardization  used  in  the 
polynomials  (A^v^(x)}  appears  to  be  quite  natural  in  many  ways.  For 
instance,  the  generating  function  defining  {A^v^(x)}  holds  for  all 
v > “1/2,  while  a separate  generating  function  for  v = 0 is  required  for 

the  polynomials  {C^(x)}.  Additionally,  the  factor  l/v  must  be  included 

to  give  the  well  known  result  that 

C°(x)  = lim  (l/v)  <£(x),  p 4 0. 
v+0  P 


Yet  with  Saw's  standardization,  no  additional  factors  are  necessary  as 


( \>)  P+v 

lim  A^V'(x)  = lim CV(x) 

v+0  p v+0  v P 


= lim CV(x)  + lim  CV(x) 

v+0  v P v+0  P 


= P 


Cp(x) 


is  seen  from 
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=A<°>(x),  p 4 0. 


For  p = 0,  both  Cq(x)  and  Aq^(x)  are  defined  to  be  unity. 

The  polynomials  |a£v^(x)}  possess  the  orthogonality  property  that 
for  p i q, 

/ a£v)(u/_§)  a(v)(u_'£)  dU(v)(_u)  = J1  a£v)(x)  A<v)(x)  dW(v)(x)  = 0. 

Sm  4 -1  P q 


Notice  that  from  the  second  integral  above  we  see  that  the  polynomials 
are  orthogonal  on  the  interval  [-1,1]  with  respect  to  the 
weight  function 


d 

dx 


W(  v\x) 


(1-x2) 


v-1/2 


B ( 1 / 2 , v+1/2) 


We  believe  the  gain  in  simplicity  and  economy  using  the  polynomials 
{Ap^(x)}  is  worth  the  slight  modification  to  the  known  formulas  for 
Gegenbauer  polynomials.  Indeed  the  greatest  gain  in  this  redefinition 
is  exhibited  by  the  formula 


Js  ApV)(u.'l)  AqV)(n_'_n)  dU(  v)(_u)  = 


0 


P 4 q , 


p = q. 


given  by  Saw  (1984).  Henceforth  we  will  use  the  polynomials  {a£v^(x)} 
exclusively,  rather  than  the  traditional  polynomials  fc^(x)}.  The  only 
dilferences  between  the  polynomials  (C^(x)}  and  {ApV^(x)J  occur  in  the 
respective  standardizations  and  the  fact  that  the  weight  functions 


35 


differ  by  a constant  factor.  Still,  the  two  sets  of  polynomials  share 
the  same  basic  orthogonality  relationship,  so  that  we  shall  also  refer 


Since  we  will  require  extensive  use  of  the  properties  of  Gegenbauer 
polynomials  in  the  sequel  and  since  it  is  cumbersome  to  collect  and 
continually  modify  certain  formulas  to  our  notation,  we  record  here  a 


For  definitions  of  other  notations  used  in  this  section,  see  the 
Appendix  listing  the  index  of  notations.  Also  for  a further  account  of 
the  results  given  in  this  section  and  for  other  results,  see  Abramowitz 
and  Stegun  (1972),  Carlson  (1977),  Durand,  Fishbane  and  Simmons  (1976), 
Erdelyi  (1953),  Gegenbauer  (1877,  1884,  1890,  1893),  Hsu  (1938),  Madhi 
and  Borah  (1984),  Rainville  (I960),  Robin  (1957),  Saw  (1984)  and  Szego 
(1975). 

2.2.1  Asymptotic  Relations 


polynomials . 


(2.2.1)  A(v)(cos9) 


2(p,v+l)  M-l  (v,m)(l-v,m) 


P 


p!  m=0  m!  (p+v-m,m) 


m!  (p+v-m,m) 


cos[(p+v-m)9  - (v4hi)tt/2] 


-M-l/2 


x 


(2sin9) V+m 


+ 0(p 


) 


0 < 9 < it,  v 4 0,-1, -2,... 
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. . p!  A^V\l)  M-l  (l-v,m) 

(2.2.2)  A(v)(cos9)  = E l 

E B(v,\H-l)  m=0  m!  (v+m,p+l) 


X 


cos  [ (p+v-hn)  0 - (v+m)ir/2] 
(2sin9) 


+ Ve) 


O<V<1,O<0<1T, 


p!  A(v)(l)  ( 1 -v , M)  A 

|Rm(0)|  < 2 "Iw" 

B(v,\H-l)  M!  (Mfv,p+1)  (2sin9) 


2sin0 


if  sin^B  > 1/2 


A = < 


COS0 


-1 


if  sin^B  < 1/2 


(2.2.3)  A(v)(cos9) 
P 


(p/2) V 1 (p+v) 
T(v+1)  sinV0 


cos[(p+v)9  - vir/2  ] 


as  p + »,  sin0  4 0 


2.2.2  Christoff el-Darboux  Formula 


. ( v)  , . ( v) , , 

n A (x)  A (y) 

(2.2.4)  2(x-y)  \ — E ?-t-E 

p-o 

P 


(tl+l)  / \ , x 

[A(^(x)  A V (y) 

( v+n+l ) A(V)(1)  "+1 

n 


-A^U)  Aj^Cy)], 


n < -2  v 
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2.2.3  Differential  Equations 


a(x)y"  + b(x)y'  + c(x)  = 0 


y 

a(x) 

b(x) 

c(x) 

(2.2.5) 

ApV) (x) 

1-x2 

-(2 v+1 )x 

p(p+2v) 

(2.2.6) 

DnA^v)(x) 

1-x2 

-(2 v+2n+l)x 

(p-n)  (p+2  v+-n) 

(2.2.7) 

(l_x2)V-i/2 

x A^v'(x) 

1-x2 

(2 v-3)x 

(p+1)  (p+2v-l) 

(2.2.8) 

sinvx  A^v^(x) 

1 

0 

v(l-v) 

(p+v)2  + 

sin  x 

(2.2.9) 

(1-X2)v/2-1/4 
x A<u,(x) 

1 

0 

(p+v)2  2+4v-4v2+x2 

1-x2  4(l-x2)2 

2.2.4  Differential  Relations 

(2.2.10)  DnA<v)(x)  = 2n  (vfl.n)  A^n)(x) 


(2.2.1  1) 


(2.2.12) 

(2.2.13) 

(2.2.14) 


DA<v)(0) 


2 (p+v)  (v+l,m)  (-l)m/m! 
0 


p = 2nrt-l 
p = 2m 


vp!(-l)P  ApV^ [x(x2-l )-1 /2]  = (p+v)(x2-l)vfP/2  dP[(x2-1)_v] 


x(p+v-l)  DA<v)(x)  - (p+v)  DA^J(x)  = p(p+v-l)  a£v)(x) 

(p+v-l)(x2-l)  Da£v)(x)  = p(p+v-1)xa£v)(x) 

~ (p+v) (p+2v-l)Ap_[ v^(x) 
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(2.2.15)  (p+v-l)DA^j(x)  - (p+vfl)DA^(x)  = 2[  (p+v)2-l  ]A<  v)  (x) 

(2.2.16)  (p+v)DAp^(x)  - (p+v+l )xDApV^ (x)  = (p+2v)  (p+v+1  )A^ v-) (x) 


(2.2.17)  (p+v+l)(l-x2)  Da£v)(x)  = (P+2v)(P+vH)xa£v)(x) 

- (p+1)  (p+v)A^(x) 

(2.2.18)  ( x- 1 ) [ ( p+v- 1 ) DAp  v } ( x)  + (p+v)A^1)(x)]  = p(p+v-l)A^  v)  (x) 

“(p+v)  (p+2v-l)Ap_[  v\x) 

(2.2.19)  (x-1)  DA^u)(x)  = pA<v)(x)  + 2(p+v)  ?l  (-l)P"k  A1(v)(x) 

P k=0  ^ 


(2.2.20)  Pa£v)(x)  + (1-x)  DA^v;(x)  = 2(p+v) 


(v), 


P-1  1 


k=0  p+k 


[^V>(K) 


+ (x-1)  DA^JU)] 


2.2.5  Explicit  Representations 

, \ (P+v)  [p/2]  (-l)k  T(p+v-k) 

(2.2.21)  A^>(x) l (2x)P_2k  , p 4 0 

r(v+l)  k=0  k!  (p-2k) ! 

/ x (x  [p/2]  (1/2, k)  xP  2k  (x2-l)k 

(2.2.22)  Ap  ' (x)  = p!  A^V)(1)  \ , p 4 0 

p k=0  (2k)!  (p-2k) ! (v+l/2,k) 

P (p+2 v, k)  [(x-l)/2]k 

(2.2.23)  Ap  ^ ( x)  = p!  A;V)(1)  l - L.  y p 4 0 

P k=0  k!  (p-k) ! (v+l/2,k) 

(2.2.24)  A<v)(x)  = p!  2~P  A(v)(l)  ? 

P k=0  k!  (p-k)!  ( v+1/2 , k)  (p+vf  1/2 , -k) 

P i 0 
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(p+v) 

(2.2.25)  A<;v)(x)  

P r(v+l) 


l 


k=0 


r(v+k)  ( v, p-k)  , /9 

[ x_(x2_1)1/2 

k!  (p-k)  ! 


x [ x+( x2- 1 ) 1 /2 ] P-k , 


p ¥ 0 


(2.2.26)  A^v)(x)  = 1 


p+v  p ( v, k)  ( v, p-k) 

(2.2.27)  A^v'(cos9)  V cos(p-2k)9  , v / 0 

v k=0  k!  (p-k) ! 

A^V\l)  00  (p+2v,k)  cos[ (p+2v+2k)  9-vir ] 

(2.2.28)  A<v)(cos9)  = — £ £ , 

B(v,v+1)  k=0  k!  (v+k,p+l) 

0<v<1,0<9<tt 


(2.2.29)  Ap^(cos9)  = 2 cos  p9, 


p = 1,2, 


(2.2.30)  A^1 ^ (cos9)  = 


(p+1)  sin  (p+l)9 


sin  9 


p 0,1,... 


r (-Dp/2  a£^2)(1) 


(2.2.31)  a£v)(0)  = 


P 0,2,4,... 
(p,v)  (0,0) 

P — 1*3,5,... 


(2.2.32)  a£v)(1)  = 


2(p+v)  r(p+2v) 
p!  f(2v+l) 


< 1 


(p,v)  t (0,0) 

p = 0 


v_ 


v = 0,  p = 1,2,... 
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Table  2.1 


Gegenbauer  Polynomials  A 


(0) 

P 


(x) 


on  S£  for  1 < p < 10 


A 

A 

A 

A 

A 

A 

A 

A 

A 

A 


<0)(*) 

= 

2x 

£0)(x) 

= 

4x2  - 2 

(°}(x) 

= 

8x3  - 6x 

£0)(x) 

= 

16x4  - 16x2  + 2 

<0)(x) 

= 

32x5  - 4 Ox3  + lOx 

= 

64x6  - 96x4  + 36x2  - 2 

$°h*) 

= 

128x7  - 224x5  + 1 1 2x3  - 14x 

(a0>U) 

= 

256x8  - 512x6  + 3 2 Ox4  - 64x2  + 2 

|0,U) 

= 

5 1 2x9  - 1 152x7  + 864x5  - 240x3  + 18x 

(0)  ( -x) 
10  w 

= 

1024x10  - 2560x8  + 2240x6  - 800x4  + 

2 


2 
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Table  2.2  Gegenbauer  Polynomials  A^^^(x)  on  S3  for  1 < p < 10 


(1/2) 


(x)  = 3x 


A^1/2)(x)  = - (3x2  _ i] 
2 


A^1/2)(x)  = - (5x3  _ 3x) 
2 


a£1/2)(x)  = - (35x4  - 3 Ox2  + 3) 
8 


A^1/2)(x)  = — (63x5  - 7 Ox 3 + 15x) 
8 


11 


(1/2), 


13 


>■6  (x)  = — (231x6  - 315x4  + 105x2  - 5) 


A(l/2) 


16 

15 


7 '(x)  = — (429x7  - 693x5  + 315x3  - 35x) 

16 

aJ1/2)(x)  = (6435x8  - 12012x8  + 6930x4  - 1260x2  + 35) 


A^1/2)(x)  = (12155x3  - 25740x7  + 18018x5  - 4620x3  + 315x) 

128  J 


128 
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a10/2)(x)  (46189x10  - 109395x8  + 90090x6  - 30030x4  + 3465x2  - 63) 

256  J 
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2.2.6  Generating  Functions 


(2.2.33) 


(2.2.34) 


(2.2.35) 


(2.2.36) 


(2.2.37) 


(2.2.38) 


l ApV>(x)  Xp  ' 
3=0  P 


l-*2 

(l-2Xx+X2) 


v+1  ’ 


N < 1 


l A<u>(x)  Ap  - (1-2AX+A2)'"  , \i  4 0,  | xl  < l 

p=0  P+V  p 


- .P  = n-^_2v 


p=0  p+v 


Ap  (x)  Xp  = (1-X)  1FQ [ v; - 


-2X(l-x) 


■]  . 

(l-X)2 

v o,  lx!  < i 


00  1 


1 + 1 A^^(x)  XP  = -ln(l-2Xx+-X2)  , 

p=l  p P 


< 1 


(v+l/2,p)  A^V\x) 


l 2 Xp  = r 1 [ (1-Xx+r)  /2  ] 


1/2-v 


p=0  p!  A<v)(l) 


f 0,  I X I < 1,  r = (l-2Xx+X2)1^2 


00  A^V\cos9)  XP 

1 -2- 


p=0  ApV\l)  p! 


T(\H-l/2)  exp(Xcos0) 

r”“r:v-i72  Jv-i/2(Xsin9)  » 


A . _ \ 

2 sm0  J 


v / 0 


00  A^V^(cos9)  XP  2 

E ■expUcos9)  0Fl[-;^l/2;~-sin29]  , 

p=0  Ap  (1)  p! 


(2.2.39) 
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00  (a,p)  A^V^(x)  a a+l  X2(x2-l) 

(2.2.40)  l 2 XP=  (1-Xx)  a F [ ;v+l/2; ] 

p=0  p!  A^v;(l)  2 2 (l-Xx)2 

00  A^V\x)  XP  X(x-l) 

(2.2.41)  £ ^ = q^i  [~>  ; ] 

p=0  ( v+l/2, p)  A^V'(l)  p!  2 

X(x+l) 

x 0F1[-;nH-1/2;~ ] 

00  ( a, p)  (2 v-a, p)  A^V\x)  t-X-r 

(2.2.42)  l XP  = „F  [ a,2v-a; v+l/2; ] 

p=0  p!  (v+l/2, p)  A',v;(l)  1 2 

P 

1+X-r 

x F [ a,2v-a;  v+l/2; ] , r = (l-2Xx+X2)l/ 

2 

2.2.7  Hypergeometric  Representations 

(2.2.43)  A<v)(x)  = A^v) (1 ) F [ -p,p+2v; v+ 1/2 ] 

» r ^ 4 2 

(2.2.44)  A^  \x)  = (-l)pA^v^(l)  _ F .[  -p,p+2v;  v+l/2 ; ] 

‘ r A 2 

, (v)  2P  (v+Up)  2 

(2.2.45)  Ap  ;(x)  (x-l)p  2Fl[  -p,l/2-p-v;  l-2p-2v; ] 

P!  l-x 

Ap  }(x)  =ApV)(l)  ( )P  2F  [ -p,l/2-p-v;  v+l/2; ] 

x+l 


(2.2.46) 


45 


(2.2.47) 

/ v (v+l,p)  (2x)P  -p  1-p  1 

;(x)  „F  [ , ; 1-p-v; ] 

P p!  2 1 2 2 x2 

(2.2.48) 

. -p  1-p  x2-l 

A(v)(x)  = A V (1 ) xp  F [ , ; vH/2 ; ] 

P P 2 11  2 2 x2 

(2.2.49) 

, x ( v+1 , p)  (x+(x2-l)  ^2  )P  . „ 

A^V  (x)  2F  J-p.v;  1-p-v;  (x-(x2-l)  ) ] 

P p! 

(2.2.50) 

A2p^(x)  = (-1)P  A^V^V)  2F  i [ "P.P+v;  1/2  ;x2  ] 

(2.2.51) 

A2p^(X)  =A2p^^  2F 1 1 -P>P+v>'r*’^/2;l-x^  ] 

(2.2.52) 

(-l)p 

A^p^i (x)  (2p+v+l)  ( vH-1 , p)  (2x)2F1  [-p,p+v+-l ; 3/2  ;x2] 

P p! 

(2.2.53) 

A2p+1^  = ^p+l^  2F i[-p»p+vfl ; v+1/2;  1_x2  ] 

2.2.8  Inequalities 


(2.2.54) 

|a£v)(x)|  < max  |a^v^(x)(  = A^(l),  v > 0,  -1  < x < 1 

P -1 <x<l  P P 

(2.2.55) 

|A2p^(x)|  < iaax  |a^(x)|  = |a2d^°^I 

-1<x<1  P P 

-1  < x < 1 , -p  < v < 0,  v not  an  integer 

(2.2.56) 

/ \ . (p+v) 

I Ap  ' (cosG)  | < (p/2) v 1 , 0 < v < 1 , 0 < 9 < tt 

r(v+l)  sinV0 
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(2.2.57)  max  |^2d+1^x^I  ^ ^/p!  [ (2p+l)  (2p+2v+l)  ] 

-1<x<1  P 

x | (2p+v+l) (v+l,  p)  | , — p— 1/2  < v < 0,  v not  an  integer 


2.2.9  Integral  Representations 


(2.2.58)  a£v\x)  = — 2 \ [ x+(x2-1)1"'/2coS(})  ]P  (sin<}>)2v  ^dtj) 

B(l/2,v)  Jn 


v > 0 


(2.2.59)  ApV^ (cos6) 


2V  A(V)(1)  (sine)1"2" 
E 

B(l/2,v) 


s 


cos(p+v)  <j> 


(cOS<j>  - cos  0 ) 


; — d. 
1-v 


v>0,0<9<7r 


2.2.10  Integrals  of  Gegenbauer  Polynomials 


b ( >.  [ (p+v-l)A^"j(x)  - (p+v+l)A^(x)]b 

(2.2.60)  / A£v)(x)dx=- Ktl 2-i -a 


2 [ (p+v)^  - l] 


(p+v)2  ¥ l 


(2.2.61)  / i ApV) (x)  dW(v)(x)  = - 


p = 0 

P = 1,2,... 


(2.2.62)  r A^Cy)  dW<">(y)  - --------  A(v|U(x)  „(u+1)U) 

_1  p(p+2u)  >r'1 


p-  1,2,...,  W^’OO W^Cx)  - a-x2)',_1/2 

dx 


/ B(l/2, v+1/2) 
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(2.2.63) 


(2.2.64) 


(2.2.65) 


(2.2.66) 


(2.2.67) 


(2.2.68) 


/ A<v)(y)  dW<v> 


,(v) 


-1 


[(■ 


1-x 


(y)  [2(v+l)w^v^(0)A^^t1^(0) 

p ( p+2  v ) 

- (2v+l)w('  v+1^  (x)A^^I1^  (x)  ] , p = 1,2,... 

) cos(0+<}))  + ( ) cos(9-<j))]  dW^V  1/2)(x) 

2 

= ApV^  (cos  9)  A^(cos<t>)  / Ap  ( 1 ) 


[ exp(iz  cos9)  A„ 

0 P 


(v) 


(cos0) 


(sin9)^v  d9  = 


tt21-v  ( v , v)  A^v^(1)  tP  z~v  Jp+V(z) 


.1 

-1 


A<v)(x)  A<v)(x)  dW(v)(x)  = 


A<V>(1) 

0 


p = q 

P ¥ q 


A^v)(a'u)A^v)(B'u)dW(v)(u)  = 

Sm  p __  q __ 


ApV)(a/l) 

0 


p = q 

P ^ q 


-q) (p+q+2v) 


2v+l 
2( v+1 ) 


fh  A^v)(x)  A^v)(x)  dW(v)(x)  = 

[w(v+1)(x)  {a£v)(x)DA<v)(x)  - A^v)(x)DA^v)(x)}]b, 
r 4 H P a 


p ¥ q 
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(2.2.69) 


(2.2.70) 


(2.2.71) 


(2.2.72) 


(2.2.73) 


I(v)(p,q,r)  = J1  a£v)(x)A<v)(x)a(v)(x)  dW(v)(x)  = 


(p+v) (q+v) (r+v)B(v,s)  / s— p+v-1 

v2  (s+v)  B(2 v, s)  \ s-p 


s— q+v-1  s~  r+v-1  \ 

s-r  ] 


s-q 


> 0,  p < q < r,  p+q  > r,  p+q+r  = 2s  is  even 


i(0)(p>q>r)  = 


p < q < r,  p+q  i r 


p<q<r,  p+q=r 


I(1)(p,q,r)  = 


(p+1)  (q+1)  (r+1) 


p < q < r,  p+q  > r 
p+q+r  is  even 

otherwise 


r 


I(v)(l,q,r)  = 


< 


r(v+i)  . 

AtV;(l) 

v+r 

(vfl)(r+2v) 
v+r 

0 


a‘u)(1) 


q = r-1 


q = r+1 


q-r  > 1 


, x p+v  q+1  ( s 

I v (p>q»r)  = [ (p-l,q+l  ,r) 

p q+V+1 


q+2v-i  p+2v~2 

+ rvMp-l,q-l,r) P vHp-2,q,r)  ], 

q+v-1  p+v-2 


p>2,  q>l,  r>0 
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2.2.10  Miscellaneous  Relations 


(2.2.74) 

a£v)(-x)  = (-l)Vv)(x) 

(2.2.75) 

, , (-l/2)Pp! (2p+2v-l) (x2-l)p/2  , . 

A(v)(x)  a‘v,(1) 

(2v-l)  (v+l/2,p) 

x a(1/2"P"v)(x(x2-1)-1/2) 

(2.2.76) 

ApV^(cos0  cos4>  + x sin9  sin<j>)  = 

P ( v+1  ,k)  (sin0  sin<j))k  ( } (v+k)  ( v-1/2) 

I t } }(cos0)  A (co s ) A (x) 

k=0  (v+l/2,k)  kK  , '(1)  p P 

p-k 

(2.2.77) 

.(l)rr  x+1  -.1/2 > 

I *«•)<,>.  taJLxU . , 

p=l  P 2n  + 1 

(2.2.78) 

n A^(x)  , . 

p=0  2p+l 

.2.11  Recurrence  Relations 


(2.2.79) 

p( p+v-2)ApV^  (x)  = 2(p+v)  (p+v-2)xA^  (x) 

- (p+v) (p+2v-2)Ap^2(x) 

(2.2.80) 

[(p+v)2-l]A^71)(x)  = u(p+v-l)A^(x)  - u(p+v+l)A^1)(x) 
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(2.2.81)  2(  v+1)  (p+v)  (1-x2)A^^^  (x)  = (p+v)  (p+2v-l )xA^^  (x) 

- p(p+v— l)ApV\x) 

(2.2.82)  2(v+l) (p+v) (l-x2)ApV+1^(x)  = (p+2 v) (p+v+1 )a£v) (x) 

- (p+1)  (p+v)xA^(x) 

(2.2.83)  (p+2)  (p+v)A^21^x)  = 2v(p+v)xA^  (x)  - 2v(p+v+l)A^(x) 


(2.2.84)  ( v-1)  [ (p+v)2-l  ]a£^  ^ (x)  = v(p+v)  (p+u-l)A^^  (x) 

- 2v[  (p+v)2-l  ]xApU^  (x)  + v(p+v)  (p+v+1  )A^^(x) 

2.2.12  Relations  with  Chebyshev  Polynomials 

(2.2.85)  Ap0^ (x)  = 2 T (x) 

Tp(x)  is  a Chebyshev  polynomial  of  the  first  kind. 

A^ (cos9) 

(2.2.86)  lim  — = T (cos9) 

v+0  A' V v 1)  P 

(2.2.87)  Ap^(x)  = (p+1)  U (x) 


up(x)  is  a Chebyshev  polynomial  of  the  second  kind. 
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2.2.13  Relations  with  Hermite  Polynomials 

, . [p/2]  (-l)k(p+v)r(p+v-k) 

(2.2.88)  A^Cx)  = l F [-k,p+v-k;-;l] H (x) 

k=0  U k!  (p-2k)  ! T(vH-l)  P"Zlc 

Hp(x)  is  a Hermite  polynomial. 

tp/2]  (-i)kr(v+i)  . , 

(2.2.89)  H (x)  = p!  I F [-k;p+v-2k+l;l] A , (x) 

k=0  k!  r(p+v-2k+l) 

*(")(x,-1/2) 

(2.2.90)  H (x)  - p!  lim  --2----- 

P p/2-1,  . x 

V>oo  V (p+v) 

2.2.14  Relations  with  Jacobi  Polynomials 

o'  A^ ( 1 ) 

(2.2.91)  a£v)(x)  = ----- e p(v-l/2,v-l/2)(x) 

P (v+l/2,p)  p 

p(a,8)(x)  a Cobi  polynomial. 

, s 2 (v,p+l)  A(vfl/2)(x) 

(2.2.92)  p(y’v)(x)  2 

P p!  (2 v+1 ) A^V+1/2) (1) 

p!  A(v/2)(1) 

(2.2.93)  A^(x)  = 2 P(v~1/2,V-1/2)(2x2-l) 

(1/2, p)  p 

(2.2.94)  A2p^(x)  = (v+1)  B(l/2,p+l)  xA^  ^ ^1)  /2  ) (i  ) p(  v-1/2  , 1/2)  ^x2_ 


1) 
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2.2.15  Relations  with  Legendre  Functions 


(2.2.95)  a£1/2)(x)  = (2p+l)  Pp(x) 

Pp(x)  is  a Legendre  polynomial. 

. s.  (p+v)  [p/2]  ( v-l/2,k) ( v,p-k) (l+2p-4k) 

(2.2.96)  A^)(x)  = l Po-2k(x)  * 

V v k=0  k!  (3/2, p-k)  P 

v f 0 

(2p+l)  2k  k!  dk 

(2.2.97)  A<£J1/2)(x) P (x)  , p > k 

P (2k+l)  (2k)!  dxk  P 

( x (p+v)  r(i/2)  r(p+2v)  . . . 

(2.2.98)  A<v)(x)  [4(x2-1)]1/4"v/2  P^2"”  (x) 

P p!  r( v+1)  P+v-1/2 

P?(x)  is  an  associated  Legendre  function  of  the  first  kind. 


2.2.16  Rodrigue's  Formula 


(-1)P  A(V)(1) 

(2.2.99)  A<">(x)  - e-TZY-  DP[(l-x2)P  w(v)(x)] 

2P  (v+l/2,p)  w (x) 


w 


(v) 


(x)  is  as  defined  in  (2.2.62) . 


2.2.17  Series  and  Expansions  Involving  Gegenbauer  Polynomials 


(2.2.100)  XP  = 


[p/2] 

r 

1 ) 

L 

k=0 


.(v)  ( N 
A (x) 

__p-2k 


k!  (v+1, p-k) 


p > 0,  -1  < x < 1 


See  Table  2.4  for  explicit  formulas  for  1 < p < 10. 
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(2.2.101) 


(2.2.102) 


(2.2.103) 


(2.2.104) 


(2.2.105) 


(2.2.106) 


(2.2.107) 


exp(xt)  = (t/2)"v  r(vfl)  l Ip+V(t)  A<v)(x) 


p=0 


v > 0,  -1  < x < 1 


sgn(x)  = l 


(-1) 


k=0  it 


B (k+1/2 , +1)  ^2k+l ’ 


v > -1/2,  -1  < x < 1 


(-1) 


k-1 


+ l 

v+1  k=  1 tt  ( 2k~  1 ) 


B(k+l/2,v+l)  ^2k^ (x)  ’ 


v > -1/2,  -1  < x < 1 


p (~l)k  A^V\x) 

(l-x)P  = 2Pp!  r(2v+l)  (v+l/2,p)  l , 

k=0  (p-k)  ! F ( p+2 v+k+1 ) 

-1  < x < 1,  -p  < (v+l)/2  if  v > 0,  -p  < v +1/2  if  -1/2  < v < 0 


-vir2v  = l 


Ilx_-Xll 


p=0  p+v  II  x II 


(v), 

p+2v  Ap  ^cos0)  » 


x,x  £ Em,  H_xll  > »X'>»  cos9  = Hxll  llx11 ) 


/ • o 

4 sin  - 


J^(2x  sin  -|  ) ® 


r(v+l)  p=0 


" l_  [Jp+V(x)]2  A^v)(cos0) 


(v)  i-xx  (v)  ” (-1)k  ^ A^(x) 

A ("f")  = P!  A V)(D  I fv) 

P k=0  k!  (p-k)!  A^V;(1) 


r = (l-2A3d-A2)1/2,  I xl  < 1 
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(2.2.108) 


(2.2.109) 


(2.2.110) 


(2.2.111) 


(2.2.112)  ( 


/ n , ,9  (p+v)  00  A^y?(x)  xk 

A(v)^  j = rP+2v___ V E+k 

P r p k=0  (p+v+k)  B(p,k+1) 

r = (l-2Ax+A2)1/2 

, . / sin9\P  p (2v+k,p-k) 

Ap  ' (cos0)  = (p+v)  I j l 

\sin<t>  / k=0  (p-k)  ! (k+v) 

/ sin(<J>-9)\  P k . 

* A^CO.*, 

\ sin0  / 

sin0  sin<j>  4 0 

/ 2p  (-l)kA^V\cos9)A^V^  (cos9) 

ApV  (cos  20)  = v(p+v)  l * ^ , 

k=0  (k+v)  (2p-k) 

( m a(v/2-1/2)(i) 

A^v/2)(cos(0-<},))  E 

v(2p+v-l ,2)  B( v/2, v/2) 

p (p-k)!  (v+2k)  ( v+2k-l  ,2)  (4  sin0  sin<|>)k 

x l 

k=0  T(p+v+k) 

(v/2+k),  . (v/2+k), 

x Ap_k  (cos0)  Ap_k  (cosij)) 

-2_)  (sin0  sin^)  ^v-l/2^a  sin4>)  exp(ia  cosB  cos<}>) 

00  i1^!  J (a)  . . 

- „2  J k+v. A,(U)(COS0)  A<v)(cos*)  , 

k=0  (k+v) (v,2k)(2v,k)  * * 


V>0,  0<9<TT,  0<(j)< 


7T 


Table  2.4  Coefficients  b for  xp  in  terras  of  A^v'(x)  on 
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2.3  Density  Representation  on  the  Hypersphere 

In  section  1.4  densities  on  Sm  of  practical  interest  were  noted  to 

be  of  the  form  f(Xu_'_£)  with  mode  and  rotational  symmetry  about 

Here  we  will  present  representations  of  such  densities  in  terms  of 

spherical  harmonic  functions  on  S . 

m 

As  a starting  point,  we  will  motivate  our  choice  of  density 
decomposition  of  Sm  by  considering  the  exponential  family  of  densities 
of  Sm.  Watson  (1983)  gives  the  exponential  family  to  be 

k 

c exP{  I Gp(Au'£)}  u.>A  e Sm 

p=l 

where  Gp(Au_'_£)  is  a homogeneous  polynomial  of  degree  p and  C is  a 
normalizing  constant  given  by 

c_l  " / exp{  l G (Au'pj  dU(v)(u). 

Sm  P=l  P 

Following  Watson  (1983),  the  constant  C will,  in  general  be  a 
complicated  function  of  the  coefficients  of  the  polynomials 
lGp(  ^H-'JI.)  Jp=l  • Hence  it  would  be  desirable  to  find  a simple 
representation  of  the  linear  space  of  homogeneous  polynomials  of  degree 
p.  Further,  we  require  the  representation  to  be  Invariant  under 
rotations  of  the  coordinate  system  which  leave  fixed.  That  is,  if 

— - where  Pis  an  orthogonal  matrix  which  leaves  _£  invariant,  then 

— A = (Pu)  = (Pu)'(Pi)  = u_'_£  and  thus  Gp(Aw_'_£)  = Gp(Au_'_£). 

Now  consider  Sp(m,_u),  _u  e Sm,  a spherical  harmonic  function  of 
degree  p in  m dimensions  as  introduced  in  section  2.1. 


It  is  well  known 
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(see  Muller  (1966))  that  there  are  (2p+m-2)  (p+m-3)  ! / [ (nr— 2) ! p ! J = 
ApV^(l),  u = m/2  - 1,  linear  independent  spherical  harmonics  of  degree  p 
in  m dimensions  which  (1)  form  a basis  for  the  linear  space  of 
homogeneous  polynomials  of  degree  p,  (2)  are  closed  with  respect  to 
rotations  of  the  coordinates  by  an  orthogonal  matrix  and  (3)  are 
orthogonal  over  Sm  in  the  sense  that 

/ S.(ra,u)  S.(m,u)  dU^v^(u)  = 0 

J S 1 ~ 

m 

for  1 < i , j < ApV^(l)  and  i ^ j . 

For  densities  restricted  to  have  rotational  symmetry  about  some 
vector  _£  e Sffi  we  have  that  if  P is  an  orthogonal  matrix  such  that  (P_u) ' £ 
= then  Sp(m,Pu_)  = Sp(m,jj).  By  standardizing  S (m,l)  to  be  unity, 
Muller  (1966)  writes  Sp(m,_u_)  as  PpCmjJl'Jj.) , a Legendre  polynomial  of 
degree  p in  m dimensions  where  PpCnijJl'j;)  is  uniquely  determined. 

On  the  circle,  S2,  Pp(2,u_'jj)  may  be  expressed  in  terms  of  the 
Chebyshev  polynomials  T (u_'J,) , or  equivalently  by  (2.2.85),  the 
Gegenbauer  polynomials  A^Cju'j;) . On  the  sphere,  S3,  Pp(3,u_'_£)  is  the 
usual  Legendre  polynomial  P (u. '_£) , which  by  (2.2.95)  may  be  expressed  in 
terms  of  the  Gegenbauer  polynomials  A^/^u^Jj) . In  general,  for  t = 
u/JL  e [-1,1],  Pp(m,t)  satisfies 

[(l-t2)D2  - (2 v+1 )tD  + p(p+2 v) ] Pp( m, t)  = 0 

which  by  (2.2.5)  shows  that  Pp(m»_u_'j^)  may  be  written  in  terms  of 

In  fact,  A^u)(t)  =a£v)(1)  Pp(2(v+l),t)  with  t e [-1,1]. 

Thus  the  Gegenbauer  polynomials  {a£  (u_’_£)  } may  serve  to  represent 
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rotationally  symmetric  densities  on  Sm  which  are  functions  of 
homogeneous  polynomials.  Since  it  is  known  that  the  Gegenbauer 
polynomials  (A^v^(x)}  are  complete  (with  respect  to  on  Sm)  for  the 

subset  of  continuous  functions  on  Sffi  which  are  of  the  form  f(u_'_§),  we 
have  the  following: 


Definition  2.3.1: 


Let  g(u.'_§.)  be  a function  on  Sm.  Then 


l i<u)  A-», 

p=0  v F 


l, 


will  be  called  the  formal  Gegenbauer  transform  of  g(u_’_§) 
where 


g(v)  . lA(v)(1)]-l 


Provided  that 


/ g(u.'l)  A^v)(u.'i.)  dU(v)(u). 

Jm 


li.  I [g(u’£>  - l 5<v)  a'u)(u'P]2  dU^Cu)  - 0, 
n-*-®  sm  p=0  P P 


then  we  will  call 


V * (v),  , 

l gp  A (u  £) 

p=0  H v 


the  Gegenbauer  transform  of  g(u_'_|), 


* 

According  to  Erdelyi  (1953),  if  g(x)  is  a continuous  function  of 
the  real  variable  x with  -l  < x < l,  then  the  Gegenbauer  transform  of 
g(x)  exists.  Pointwise  convergence  of  \ ^ g^v^  ApV^(u_'jj)  to  g(u'£) 
almost  everywhere  with  respect  to  requires 


CD 


l 

p=0 


[gpV>]2  ApV)(t) 
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We  now  present  a table  of  the  coefficients  gpV^  for  several 
well  known  densities  on  Sm  (refer  to  Saw  (1984)).  For  the  Gegenbauer 
transform  of  other  functions  see  the  listing  of  formulas  in  subsection 


2.2.17. 
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Table  2.5  Coefficients,  gpV\  of  Gegenbauer  Transforms 


of  Densities,  g(Au_'_§),  on  S 


m 


Distribution 

g(Au'!) 

g£V\  p=  0,1,2,... 

Langevin 

Ay-l(A)  exp(Au_'_5) 
AV(A)  = qF1  (-;  v+1 ; X2 / 4 ) 

Ip+u(A)/Iv(X) 

Saw' s 

2 00  (2Au' £)k( v+l,k/2) 

i 

k=0  k! 

APM(  p/2 , p+ v+1 ; - A2 ) 
(v+l+p/2,p/2) 

Saw' s 

Hyperspherical  t 

K(V)(X) 

s 

AP(  v+s  ,p)  K^v)(A) 

(l-2Au,^+A2)Vfs 

(v+l,p)  K^U) 

v)  (A)  ]_1=2F  l (s  , v+s;  v+1 ; A2) 

Cauchy 

1-A2 

(l-2Au,^+A2)Vfl 

aP 

Cardioid 

1 + AA1(v)(u_'_^) 

l if  p = 0 
A if  p = 1 
0 if  p > 2 

Wrapped  Normal 

i »p2 

p=0  p 

AP2 

Point 

I ApV)(u'iP 

p=0  p 

1 

CHAPTER  3 

THE  STATISTIC  Yp(») 

3.1  Introduction 

In  section  1.6  tests  of  uniformity  proposed  by  Ajne  (1968)  were 
discussed.  Ajne's  tests  were  developed  for  testing  uniformity  on  the 
unit  circle  S2»  Of  central  importance  in  those  tests  is  the  statistic 
N(tt,x)  which  counts  the  number  of  sample  points  in  some  semicircle 
corresponding  to  the  interval  (x,x+irj. 

The  analogue  of  the  statistic  N(tt,x)  on  m-dimensional  spheres 
should  count  the  number  of  sample  points  on  some  m-dimensional 
hemisphere.  We  denote  the  m— dimensional  hemisphere  H ( u ) , rotationally 
symmetric  about  u_  e Sm,  by 


= l*  e Sm  : — > °}- 


Then  given  the  random  sample  x^ , . . . ,xn  from  Sm,  define  N(u_)  to  be  the 
number  of  sample  points  in  Hm(_u) . In  terms  of  indicator  functions, 

N<a>  ‘ j,IfVH>l(SJ>' 

Under  the  hypothesis  of  uniformity,  N(u)  follows  a binomial  distribution 
with  parameters  n and  p = 1/2. 

Generalizations  of  Ajne's  results  to  m-dimensional  spheres  have 
been  recounted  in  section  1.6.  However,  the  test  statistics  in  those 
generalizations  have  complicated  asymptotic  distributions.  In  section 
3.3  we  shall  propose  a test  statistic,  Y (•),  motivated  by  N(u),  which 


has  a normal  asymptotic  distribution  under  the  hypothesis  of  uniformity. 
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3.2  The  Functions  X(«)  and  6(») 

Under  the  null  hypothesis  of  uniformity  we  standardize  N(u),  the 
number  of  sample  points  _x^  » • • • which  occur  in  Hm(_u ) , and  define 

XU)  = 2n~1/2  (N(u)  - n/2). 

Both  Ajne's  statistic  I and  Prentice's  m— dimensional  analogue  A 

° m,n  ’ 

discussed  in  section  1.6,  are  functions  of  X(_u) . 

In  future  work  we  will  find  alternate  expressions  for  X(u) 
useful.  As  a first  step  toward  deriving  alternate  expressions  for  X(u), 
notice  that  to  count  the  number  of  sample  points  in  Hm(u)  is  equivalent 
to  subtracting  the  number  of  sample  points  not  in  Hm(_u)  from  the  total 
sample  size  n.  Supposing  that  no  sample  points  fall  on  the  boundary 
between  H^U)  and  Sm  - HmU)  = Hm(-u) , then  those  sample  points  not  in 
belong  to  Hm(—  u ) and  hence  N(jj_)  = n - N (— u ) . Thus  we  may  write 

X(_u)  = n ^2  ( 2N(_u_)  - n) 

= n_1/2  [N(u)  - (n  - N (u) ) ] 

= n~l/2  (N(_u)  - N (— u)  ) . 


Observe  that  N(_u)  N(  _u)  counts  the  difference  between  those  sample 
points  in  H^u)  and  those  in 

Uow  for  u_  e Sm  and  _£  e Sq  define  the  signum  function  sgn(u'f;)  = 
6(uU)  as 
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(3.2.1) 


<5(u/j,)  = 0 


-l 


1 


if  _u > 0, 
if  u_'_£  = 0, 
if  u_'_£  < 0. 


Thus  6(u_'_x)  = l if  x_  e Hffl(up  and  <S(u.'20  = “l  if  .x  e Hm(-u_).  We  may  now 
write 


J“* 

Since  N(^i)  ~ Binomial  (n,l/2)  under  the  null  hypothesis  of 
uniformity,  it  is  obvious  that 

EX(u)  = 0 and  Var(X(u))  = l. 

To  compute  Cov(X(u^ ,X(w) ) we  first  prove  the  following: 

Theorem  3.2.1:  Let  x U(S_)  and  u.w  e S . Then 

———————  — m — — — m 

E6(x_'_u)  = 0, 

Cov(  , 6(_x' w)  ) = (2/ it)  sin-l(u/w). 

Proof : Since  _x_  ~ U(Sm)  we  have 

E6(x.'_u)  = P(x_  e Hm(u.))  “ P(x  e H^-u)) 

= 1/2  - 1/2 


n 


N(u)  - N(-jj)  = ][  6(u'x. ) 

j=l  " "J 


so  that 


(3.2.2) 


X(u)  = n 


0. 
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To  compute  the  covariance  we  first  divide  Sm  into  the  following 
four  regions: 

R1  = A Hm(w)  = [x  e sm  : 2L'iL  > 0.  > °}» 

r2  = A ^(-w)  = {x_  G Sm  : _x'_ u > 0,  _x'_w  < 0 } , 

R3  = A H^Cw)  = [x  e Sm  : x.'u_  < 0,  x'w  > 0} , 

r4  = A = {*.  e sm  : x?u_  < 0,  jc'_w  < 0}. 

These  four  regions  are  illustrated  in  Figure  3.1  for  the  case  m = 2. 


Figure  3.1  Regions  Rp  R2,  R3  and  R^  for  m = 2 
Let  = P(x.  e R^)  . i=  1,  . . . ,4,  x_  ^ U(Sm)  . It  is  sufficient  to 
only  compute  because  for  x ~ U(Sm),  pl  = p4  and  p2  = P3  = 1/2  - pp 
Then  using  Theorem  1.3.1  to  write  x U(Sm)  as  x = z / !l z S where  z ■» 
Nm(0,Im)  we  have 

Pl  = p(i  £ Ri^ 

= P(_x'_u  > 0,  _x'_ w > 0) 

z ' u z ' w 

= p( > 0, > 0) 


!!  z II 
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= P(z'u  > 0,  z’w  > 0) 


= P(y1  > 0,  y2  > 0) 


where  y^  = _z_'u_  and  y2  = _z_' w_.  But 


m 


y = z'u  - N(0,  l u = l)  , 
1 i=l 


and 


ra 


y = z'w  - N(0,  l w = l)  . 

i=l 


Further, 


Cov(y1,y2)  = Cov(_z_'_u,  _z_'_w) 


m 


m 


= covf  y u.z . , y w .z . 

1-1  11  j-1  J J 


m 


= y u. w.  Cov(z . , z . ) 

.,ii  l i 

i=l 


= u'w. 


Thus , 


where  p = u'w  gives 


^ r°°  r00  ^ 2 2 

P = I V exp| (y  _ 2 py  y + y2)}  dy  dy 

2-rr  (1-p2 ) 1 J0J0  2(l-p2)  1 1 1 L V 1 
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Now  let  y^  = s(l-p2)l^2  + pt  and  y2  = t.  Since  the  Jacobian  of  the 
transformation  from  {y1>y2}  to  { s , t } is  J{ypy2  s,t}  = (l-p2)1//2,  then 


l C °°  C 00  ~i 

, = \ \ exp{ (s2  + t2)}  dsdt. 

2tt  J 0 J -p  2 


(1-P2) 


1/2 


2 2 2 

Changing  to  polar  coordinates  with  r = s + t , s = r cos6, 
t = r sin6  and  j{s,t  -►  r , 9 } = r,  we  have 


1 pn/2  r 

2tt  J -1  p J i 

tan  77~,7ui 


r exp{-r2/2}  drd9 


(1-P2) 


1 r ' -1 

| tan  -r7r 

2tt  2 (1-P2)  7 


1 1 -i  “P 

— + tan  

4 2tt  (l-p2) 


1/2 


1 1 TT 

+ { cos  Ap  } 

4 2tt  2 


-1 


1 1 -1,  . V 

— - cos  (u'w) . 

2 2ir 


Observe  that. 


S(x'u)  6(x'w)  = « 


1 if  _x  e Rj  orjc  e L, 


-1  if  x e R„  or  x e R„ 


_ “2  or  _x_  e < 
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Hence, 

Cov(6(x_'u_)  6(2c_'_w)  ) = E ( 6(_x ’u_)  5(x_'w)  ) 

= p(x^  e Rj  or  x.  e R4)  “ P(x.  e R£  or  e Rg) 

= Pt  + P4  - P2  " P3 
= 4Pl  - 1 

1 1 -1 

= 4 { cos  1 (u'w) } - 1 

2 2 TT 

2 

= 1 - COS  (u'w) 

IT 

2 , " 

= l { sin  (u'w)} 

it  2 

2 

= sin  (u'w). 

n 


With  Theorem  3.2.1  we  may  prove  the  following: 


Theorem  3.2,2:  Let  _Xj 

n-l/2(2N(u_)  - n), 


- i.i.d.  U(Sm),  j = l,...,n  and  X(u)  = 


Then 


EX(u_)  = 0 
2 

Cov(X(u)  ,X(_w) ) = - sin  l(u.'w). 

IT 

Proof:  It  has  already  been  established  that  EX(_u)  = 0.  To  compute  the 
covariance,  we  express  X(_u)  in  terms  of  the  signum  function  and  use 


Theorem  3.2.1  to  find 
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n 


Cov(X(„),XCw»  l «8'5i>'  «s’5,» 

1=1  J = 1 


n 


n 1 l l E(5(u'x. ) 5(w'x  )) 
i=l  j=l  ‘ J 


n n 


n 1 l E(6(u_,3ti)  SOw'iLp) 
i=l 


+ n"  ll  E(5(_u’xi))  E ( 6 (w ’ x.j  ) ) 


n 2 _ 

= n 1 l sin  vu/w) 

i=l  it 


sin  *(u'w). 


Returning  to  the  signum  function  S(u_'_£),  we  verify  formula 
(2.2.102)  that  the  Gegenbauer  transform  of  6(u'5)  is 


l «. 

p=l 


where 


(_! ) (p— 1) /2 


B (p/2 , v+l ) 


6p" 


P 1,3,5,..., 


p 0,2, A,.... 


By  definition  2.3.1  and  formula  (2.2.67), 
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/ 5(u'0  A<v)(u'0  dU(v)(u) 

sm  — p — ~ 

in 


- / ( l & A<v>(u'{)}  a£u)(u'{)  dU(u)(u) 

S„  q-0  ’ q P 


- 1 / (A<“>(u'5)l2  dD^Cu) 

P ' C 1 P 

m 


■ «p  A<v)Cl). 


Since  A^v^(x)  is  an  odd  (even)  function  for  p odd  (even)  and  6(x)  i 
odd  function,  we  have 


s an 


«P  Ap',)<1>  ■ < 


2 / A^v\x)  dW(v)(x) 

0 p 


P 1,3,5, ..., 


p = 0,2,4, ... . 


Then  for  p = 1,3,5,...,  formula  (2.2.63)  gives 


<5„  = 


4(  v+1 ) A(y|1}(0) 

.2ZL 


p(p+2v)  B ( 1/2, v+1/2)  A (v(l) 

P 


Upon  simplification  using  formulas  (2.2.31)  and  (2.2.32)  we  get 


(_1}(p-l)/2 


6P  = 


B(p/2,  vfl)  , 


P = 1,3,5, 


We  now  turn  to  the  Gegenbauer  expansion  of  Cov(X(_£)  ,X(_n)  ) = 

2/tr  sin  (_§_'jt)  = x(w),  say,  where  w = e Sm*  We  note  in 

1/2 

passing,  that  x(w  ) is  the  distribution  function  of  the  well  known 


arcsine  distribution. 
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Theorem  3 .2.3;  The  coefficients  of  the  Gegenbauer  transform  of  x(w)  = 
2/ tt  sin  ^w  are  given  by 


'p-3 


P = 0,1,2, 


where  6p,  p = 0,1,2,...  are  the  coefficients  of  the  Gegenbauer 


transform  of  the  signum  function  6(w) . 


Proof : Let  w = with  e Sm.  Then 


t(w)  = x(_£’_n) 


= Cov(6(u,C),5(u'n)), 


“ ' 


- El(  l_  «p  A<v)(u\£))(  l 6a 


p=0 


q=0 


q q 


= E{  I.  l_  5q  A^v)(u-0  A^v\u_'_n)} 


p=0  q=0 


p=0 


by  (2.2.67) 


' l.  V ApU,<i.’jl)- 


p=0 


Thus  by  the  uniqueness  of  the  coefficients,  we  have  xp  = 5p, 
p 0,1,2,.... 
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3.3  The  Statistic  Yp(*)  under  the  Hypothesis  of  Uniformity 

In  this  section  we  shall  introduce  the  statistic  Yp(»)  and  note 
several  properties  of  Y (•)  under  the  hypothesis  of  uniformity.  Section 
3.4  will  examine  Yp(»)  under  various  alternative  hypotheses  while 
Chapter  4 will  discuss  testing  uniformity  with  Yp(»). 

The  generalization,  A , of  Ajne’s  statistic  1/ (2irn)  was  given  by 

HI  y II 

formula  (1.6.3)  and  may  be  written  in  terms  of  the  measure  as 

A , -1  / { N(_u_)  - n/2}2  dU(v)(u) 

m*n  Sm 

= 1/4  / (X(u_)}2  dU(v)(_u). 

^m 

However,  it  is  somewhat  difficult  to  work  with  the  asymptotic  null 
distribution  of  Am  n which  involves  an  infinite  series  of  weighted  chi- 
square  random  variables.  We  are  thereby  motivated  to  find  a kernel  for 
the  integral  of  (N(_£)  - n/2}  which  simplifies  the  asymptotic 
distribution.  The  relation  between  this  new  integral  and  A is 
indicated  in  section  4.1.  Now  we  make  the  following: 

Definition  3.3.1:  Let  _x^ , . . . .^  be  a random  sample  from  Sffl.  Set  X(_u)  = 

n l//2(2N(u_)  - n)  where  ji  e Sm,  N(u_)  = J\(l  + 6(u'x_j))/2  and  6(w)  is 

the  signum  function.  Then  for  p = 0,1,2,...  and  e Sffl,  define 

Vi)  = (-l)(P_l)/2  J X(u)  A^Cu’O  dU(v)(u). 

m 

A computational  form  for  Yp(_^)  will  now  be  derived  by  using  the 

expression  (3.2.2)  for  X(u_)  written  in  terms  of  the  signum  function.  In 
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addition,  the  computational  form  of  Yp(_£)  will  show  that 

Yp(_£)  = 0 for  p = 0,2,4,.... 

Substituting  X(u_)  = n 6C_u_*2£j  ) and  the  Gegenbauer  transform 

of  6(w)  into  the  definition  of  Yp(_£)  yields 


Y„(£)  - (-l)(p_l)/2/  n_l/2  l 5<u_'x . ) A*v,(_u'_5_)  dU(v)(u_) 


n «>  _ 


(_U(p-l)/2  n-l/2  11(1  «»<»><«•*)) 

J-l  Sm  q=0  J P 


dU(v)(u) 


(-DCp-D/2  n-l/2  7 £ j A(v)(u'xi)A^v)(u*l)  dU(v)(u) 

j-l  Sm 


P j-l  P J 


By  expanding  6 


P 


we  have 


(3.3.1)  Yp(jj)  = 


.-1  „“l/2 


n"i/z  B(p/2,vfl)  l A<v)(£'x,) 

j-l 


P - 1,3,5, 


P = 0,2,4, 


Under  the  hypothesis  of  uniformity  we  now  examine  the  mean  and 
covariance  for  Yp(_£)  in  the  following: 
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Theorem  3.3.  1 : Let  _Xj  be  i.i.d.  U(Sm) , j - 1,  . . . ,n  and  £ Sffl. 


Then 


EYU)  - 0 


p — 0»1»2,..«, 


Cov(Yp( O ,Yq(n)  = ' 


tp  A^Cl'n) 


p = q = 1,3,5, 


otherwise, 


where  ft  } are  the  coefficients  in  the  Gegenbauer  transform  of 

1 Pp=0 

t(w)  = 2/tt  sin  ^w. 


Proof:  Since  Yp(_£)  = 0 for  p even,  we  need  only  consider  p odd.  By 

formula  (2.2.61)  and  the  computational  form  of  Y (_£)  given  in  formula 
(3.3.1)  we  get  for  p odd 


EY  (£)  = tt"1  n_1/2  B(p/2 , v+l)  l EA^^U'x.) 

j = l 


= 0. 


For  the  covariance,  we  express  Yp(jj_)  and  Yp(_n)  in  computational 
form  by  formula  (3.3.1)  and  then  apply  formulas  (2.2.67)  and  (2.2.61), 
Thus  for  p and  q both  odd, 


Cov(Yp(0,Yq(ji))  = E(Yp(J_)  YqU)) 


Eta 

2 


- 1 


- (-1)  2 V1  6 5 l l ECA^Cl'iPA^Wxj)) 


P "l-l  J-l 
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= (-D 


2 +a 

2 


- 1 


-1 


V,(  j.  E(A‘V)(i,*1)A<V)(n'xi)) 


i=l 


+ l l E(a£v)(£'x,))  E(A<v)(_n'x .))} 

i*j  4 J 


(i'jl) 


■ «p 


ApV>(i,Jl) 


where  tp  = 6p  by  Theorem  3.2.3. 

It  has  already  been  shown  that  for  p = 1,3,5,...,  Yp(jj)  = y^ 

where  y..  = (-1)(P“1)/2  n-l/2  j = l n.  Ifcreover,  from 

the  mean  and  variance  of  Yp(j[),  p = 1,3,5,...,  under  uniformity  we  see 
that  Eyj  = 0 and  VarCy^)  = Tp  n_l  a£v)(1),  j = l,...,n.  Therefore  by 
the  central  limit  theorem 


( 


T 

P 


A^U)  ),/2 

O ' 


n 


-i 

n 


A^U)  )l/2 

n J 


converges  to  N(0,l)  under  the  hypothesis  of  uniformity.  Consequently, 
we  have  shown  the  following: 
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Theorem  3 .3.2:  Under  the  hypothesis  of  uniformity,  the  asymptotic 
distribution  of  Yp(_£),  p = 1,3,5,...,  is  N(0,fp  Apv^(l)). 

The  exact  distribution  of  Yp(_£)  is  difficult  to  find  except  in  the 
special  case  m = 3 and  then  only  for  p = 1.  In  particular,  we  have  for 
P ■ 1. 


Vx> = n_1/2  6i  l Aiv)av 

j=i 

= 2 (vfl)  n_1/2  l 

j = l 

=2  (v+1)  n~1/2  &l  y 


where  y = V . _£'x. , j = l,...,n.  To  find  the  distribution  of  Yi(jj)  it 
is  sufficient  to  know  the  distribution  of  y.  As  in  section  1,5  we  may 

define  _s_=  J . Xj,  j = rn  = l[sjl  and  Xq  = _s_  / rn.  Hence 

y = 1'  I j iLj  = JL'±  = rn  i.'x 0,  j = 1,  ...,n. 

Under  uniformity  it  is  well  known  (see  Watson  (1983)  for  instance) 
that  Xq  - U(Sm)  independently  of  rn«  Put  w = so  that  by  Theorem 

1.4.1,  w - ^ j ( v+1/2 , v+1/2) . Thus  we  seek  forms  for  the  density  of 

y = rn  w. 

To  have  a one-to-one  transformation  we  define  z = w and  note  that 
the  Jacobian  of  the  transformation  from  {r,w}  to  {y,z}  is  given  by 
j{r,w  * y,z}  = z *.  We  may  now  express  g(y) , the  density  of  y,  as 

g(y)  = / h (y/z)  f ( z ) z_1  dz 
y/n 

or  equivalently  with  u = y/z, 


g(y)  - / h (u)  f ( y/u)  u 1 du 

y 
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For  a random  sample  from  a uniform  distribution  on  Sm,  hn(r) » the 
density  of  the  resultant  length,  is  known  to  be 


hn(r)  = tV+1(2vr(v+l))n_1/*9v(l“n)+1(Jv(0))nJv(9r)d0 


for  0 < r < n.  Likewise  for  a random  sample  of  observations  uniformly 
distributed  on  Sm  the  density  of  w = is 


f (w) 


(1-w2) 


v-1/2 


B ( 1 / 2 , v+1 /2) 


for  -1  < w < 1. 

For  m = 3,  jj_  e and_Xj  i.i.d.  U(Sm) , j = l,...,n,  we  have  by 
Theorem  1.4.1  that  _§_'x_j , j = l,...,n,  are  i.i.d.  uniform  [-1,1].  Thus 
y,  the  sum  of  n i.i.d.  uniform  [-1,1]  random  variables,  has  density 


g(y)  = 


1 [ (y+n) /2  ] 

l (.1}J 

2 (n-1)!  j=0 


(y+n-2j)n  1 , 


-n  < y < n and  [x]  is  the  greatest  integer  function 
the  density  for  Yj(_£)  = m n-1//2  6!  y to  be 

c r (ct+n)  /2  1 

f (t)  l (_1} 

Y1<D  2nu-l)!  jio 


Using  g(y)  we  find 


t+n-2j)n  1 


y 


for  -n/c  < t < n/c  where  c = { (nir) 1 /2r(  v+3/2) } / |2r(v+2)}. 

If  the  modal  vector  _£  e Sm  is  unknown,  we  may  estimate  _£  by  _Xq. 
Consequently,  the  distribution  of  is  also  of  interest.  Again  the 

distribution  is  difficult  except  for  p = 1 where 
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yi(xq)  = n_1/2  «i  l A J (_Xq  'x.j  ) 

j=l 


n 


(v). 


= m n 


-1/2 


n _ 


«l  ^ JEo'ij 


a n 


■ m n'1/2  si  (r/  j.  2i yi  ^ 

i=l  J=1  J 


n n 


= m n_1/2  5j  r;1  £ I ^ x 

i-1  j-1  J 


= m n"1/2  5j_  r"1  r2 


= m 


„-!«  «1  rn* 


Hence  the  density  of  Y^(xq)  may  now  be  found  from  the  density  hn(r) 


to  be 


f (t)  = c(ct)vfl[2v  r(vH)]n_1 

V*0} 

X /”9v(1-n)+1[Jv)(e)]nJv(0ct)  d9 

for  t > 0 and  where,  as  before,  c = { (m;)  ^2r(  v+3/2)  } / (2r(\H-2)}. 

O 

The  asymptotic  null  distribution  of  Y|(xq)  may  be  found  from  the 

o j o 

well  known  result  that  under  uniformity,  m/n  r^  ^ xm  as  n -*■  With 
m = 2(v  + 1)  we  have  Y2(xq)  = m2/n  52  r2  ^ m 62  as  n ^ 

Before  proceeding  with  examining  Yp(_£)  under  the  alternative 

9 

hypothesis  of  non-uniformity,  we  cite  a relation  between  Y^(jj)  and 
Y2(xq).  We  have 
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tn  / Y^(£)  dU(v)(£) 

Sm 


A ^ v)  (_i.'2Li  )A  i V ) (-5.'2Lj ) 


du(v)(i) 


n 


l 

j = l 


A\ v)  (-2-1— j ) ) 


m 0 n n 

- 67  l l x!x. 

n 1 i=l  J-l  "1_J 


m -2  2 

" S1  rn 
n 


- ^(Xq). 


3.4  The  Statistic  Y (•)  under  the  Hypothesis  of  Non-Uniformity 

Before  turning  to  testing  uniformity  in  Chapter  4 with  the 
statistic  Yp(_£),  we  consider  some  properties  of  Yp(_£)  given  the 
alternative  hypothesis  that  the  true  underlying  density  is  non- 
uniform.  In  particular,  we  take  the  alternative  hypothesis,  Hp  to  be 
that  the  density  f(_£'.x)  on  sm  is  non-uniform,  rotationally  symmetric 
about  _£  e Sm,  and  an  increasing  function  of  Jj/x_,  2L  e Moreover,  let 

l £p4u)(i'i) 

P-0 


be  the  Gegenbauer  transform  of  f(E'u). 


79 


Theorem  3.4.1;  Suppose  _Xj ,3^  are  i.i.d.  random  vectors  on  Sm  with 
common  non-uniform  density  f(_£'_x)  as  described  above.  Let 
be  the  signum  function  and  let  fp  and  <$p,  p = 0,1,2,...,  be 
the  coefficients  of  the  Gegenbauer  transforms  corresponding  to 
f(_£'x_)  and  6(_£'x_).  Then 

Y ( O = (-1)(P-1)/2  n-1/2  5 l 
V j=0  P J 

has  mean  and  covariance  given  by 

Pp  - (-1)(P-1)/2  n1/2  5p  fp  A<v)(l) 

and 

. (p+q)/2  _ , . 

Ppq  ■ <-D  * 6_  { I C r 'Cp,q,2r) 

v 4 r=0 

- *p  f,  ApV)<1)  AqU,(1)l 

« 

respectively,  where  p,q  = 0,1,2,...  and 

I(v)(p,q,r)  = / 1 A<v)(x)A(y)(x)A^^(x)  dW<v>(x) 

is  found  from  formula  (2.2.69). 

Pr oof : Taking  expectation  with  respect  to  the  non-uniform  density 

f (_?.T  JLj  ) » j“l>  •••,!!,  we  find 


EYpU)  = (-1)  (P  1)72  n“1/2  6p  l Ea(v)(_£'xj) 


= (-D<P  D/2  n-l/2  6p  ^ / A^v)U'x_.) 


j-1  m 


l fq  Aq  V)  (JL'iL-j  ^ dU(v)(x_,) 
q=0  4 4 J J 


X 
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= (-i)(P"l)/2  n-l/2  ^ j fp  ApV^ (1 ) 
= (-i)(p-D/2  nl/2  ^ fp  A£v)(l). 


To  compute  the  covariance,  apq,  we  first  consider  E(Yp(_£)  Y (_£)) 
where  the  expectation  is  again  with  respect  to  the  non-uniform  density 

j = l,...,n.  Since  Yp(jj)  = 0 for  p = 0,2,4,...,  then  obviously 
opq  = 0 if  either  p or  q is  even.  Thus  for  p,q  = 1,3,5,..., 


E<yi.)  yi))  - E{  (-D 


Eta 

2 


- 1 


n'1  «p  «.  I l 

F 4 i-1  j=l  4 J 


n n 


Eta  -i 

- (-1)  2 n_i  5p  6q  E{  ^ 


1*1  J 


- (-1) 


Eta  _i 

2 -1  T T 


*P  *q  { l fS  ApV)(i,ij)AqV>^,ij) 
j — 1 m 


00 

x l f a£v)U'x_ i)dU(v)(x_.) 
r=0  J J 


+ l l E(A Jv)(i.'Xi))  E(A<v)(£'xi))} 
i*j  4 


As  in  the  derivation  of  Up  above,  E(A^v)(_£'xi) ) = f a£u)(1).  Thus 
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E(yi>  yj.))  - (-D 


Eta 

2 


- 1 


a"1  6 6 

p q 


(Ilf, 

j-1  r=0 


T(v) , . 

I (p,q,r) 


y i 


f f A<v)(l)  A^V)(1)  }. 

p q p q J 


Using  formula  (2.2.69)  and  the  fact  that  p and  q are  both  odd,  r must  be 
even  and,  further,  r < p+q.  Hence 


E<yi>  vi»  = (-D 


Eta 

2 


- i 


n'1  6 6 

P q 


(p+q) /2 

(n  l 

r=0 


f2r  I(v)(p,q,2r) 


+ n(n-l)  fp  fq  a£v)(1)  A<v)(1)} 


= (-1) 


Eta 

2 


1 (p+q)/2  , 

« « { l f2r  rW(p,q,2r) 

H 4 r=0 


+ (n-1)  fp  fq  A<v)(l)  A<v)(l)}. 


Then  using  the  formula  for 


Up  given  above, 


°pq  - E<yi>  Vl»  - un  U 


p Mq 


= (-D 


Eta 

2 


- 1 _ _ (p+q)/2 

{P  6q  1 l y lVV'<P>q>2r) 

H 4 r=0 


r(v), 


+ (n-l)  fp  fq  A<u)(l)  A<u>U) 


n *p  ApU)d>  A<v)(l)} 
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= (-D 


Eta 

2 


- 1 


6 6 { 


p q 


(p+q) /2 

l 

r=0 


'2r 


I(u)(p,q.2r) 


ApU,(‘)  A$V><»}- 


For  the  variance  of  Yp(_£)  we  use  the  notation 
3.4.1  that 


o 

Op  and  find  from  Theorem 


a 


2 

P 


6P  ( i f2r  l(V)(P>P>2r>  - 
r=0 


Apv)(‘>  ]2 


1- 


Finally  we  consider  the  asymptotic  distribution  of  Yp(_£)  under  non- 
uniformity in  the  following: 

Theorem  3.4,2:  Let  f ( ^ *x ) be  a density  on  Sm  which  is  non-uniform, 

rotationally  symmetric  about  e Sm,  and  an  increasing  function  of 
JL'iL*  Suppose  jcp  ...  ,3^  are  i.i.d.  random  vectors  on  Sm  with  common 
density  f(£'x).  Then  Y (£)  = (-l)(P_l)/2  -1/2  y A(v)(f'x.) 
is  asymptotically  distributed  N(pp,cjp)  where  Pp  and  are  given  in 
Theorem  3.4.1. 


Pr oof . As  in  the  derivation  of  Theorem  3.3.2,  we  may  write  Y ( ^ ) = V.y. 
where  y.  = (-i)(p-D/2  ^ n~l/2  A^^’xj),  j = l,...,n.  Similar  to  the 
calculations  in  Theorem  3.4.1,  we  take  expectations  with  respect  to 

to  find  Eyj  = pp/n  and  Var(yj)  = Op/n.  Then  by  the  central 
limit  theorem 


YU)  ~ K [ ( y.  - p n"ll 
P P ^ W3  P J 


1/2  , -1/2, 
n (°p  n ) 
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converges  in  distribution  to  N(0,1).  Thus  we  deduce  that  Yp(_£)  ^ 

N(up,Op)  under  non-uniformity. 


CHAPTER  4 

TESTING  UNIFORMITY  WITH  THE  STATISTIC  Yp(  • ) 

4.1  Introduction 

In  section  3.3.  the  statistic  Y (jj)  = (-1)  ^ n 5p 

x £ . 3 = !•>•••  >n>  was  introduced.  Section  3.3  and  3.4 

J P J 

found  the  asymptotic  distribution  of  Yp(_£)  to  be  normal  under  uniformity 
and  for  non-uniform  distributions  fO^'jt)  which  are  rotationally 
symmetric  about  jj_  e Sm  and  increasing  functions  of 

Here  we  will  discuss  using  the  statistic  Yp(_£)  to  test  for 
uniformity  on  Sm.  Indeed,  we  will  see  how  Ajne's  statistic,  generalized 
to  Sm,  is  a function  of  {Yp(_£)}p>0. 

A generalization  of  Ajne's  statistic  I/(2irn)  given  by  Saw  (1984) 
was  recorded  in  formula  (1.6.5)  as 

=„(u)  -Li  I [♦(sJi)  - Ml2  ™M<-V  • 

b A=]  J 

m J 1 

Set  ) = 1 + 6 (_§_  *2Lj ) » j = where  6(w)  is  the  signum 

function.  Ihen 

n~1(4V)  = n~1L  { I 5(x’i)}  2 dU(v)(£) 

Sm  j-1  J 

= f {x(i)}2  dU(v)U) 

sm 

= E_ ({X(i.)}2).  i.  £ U(Sm) 
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where  XU)  = n V . fi(£'x.)  is  as  given  in  formula  (3.2.2).  Also 

J J 

since  X(_£)  = 2n  ^2(NU)  - n/2),  we  see  Chat  = 4n  Am  n in  which 

Affl  n is  the  generalization  of  Ajne's  statistic  to  Sm  given  by  Prentice 
(1978).  Note  that 


XU)  = n~1/2  l SU'xJ 

j-1 

= »*‘/2  l l 5pa‘v)(£*x) 

J-1  P=0  p p J 

= l (-l)(p-1)/2  Y (£)  . 

p=0  p 


Thus 


2+3  _! 

{*(£)}  " l {Y  U)}2  + l l (-D  2 YU)  YU). 


p=0 


p*q 


But  by  formula  (2.2.67),  E?(YpU)  YqU))  = 0 for  p 4 q and  - U(Sm). 
Therefore  with_£_  *■  U(Sm), 


Gn(U>  - " 

= n l E.({Tf  (S) }2) 

P=0  a p 

= n I L (yj.)}2  dU(v>(u) 

P=0  Sm  P 

so  that  Gn^v^  and  {Y2(_£)}p>g  are  functionally  related. 
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4.2  Optimal  Properties  of  Yp(») 

Optimal  properties  of  the  statistics  Y^(jj_)  and  Y^Cx^)  for  testing 
uniformity  on  Sm  arise  from  the  relations  described  in  section  3.3, 
namely , 


Y^i.)  = m n 1/2  6j  rn 

and 

^(xq)  = m n-1//2  rQ. 

Ihus  Y^(xq)  enjoys  the  optimal  properties  of  the  Rayleigh  test,  while 
Y^(j[)  shares  the  optimal  properties  associated  with  rn  'xq  = rn  cos<j> 
where  <p  is  the  angle  between  jj_  and  _Xq . Although  these  optimal 
properties  are  widely  known  for  the  circular  and  spherical  cases,  we 
shall  present  the  results  for  general  m- spheres , m > 2,  whenever 
possible. 


4*  2. 1 Uniformity  versus  Langevin  Distribution  for  Mean  Direction  Given 

Consider  the  observation  u,,...,u  e Sm  which  are  i.i.d.  with 
common  Langevin  distribution  R(v)(X,_£)  where  is  known.  We  propose 
testing  Hq : X = 0 against  H^:  X > 0 in  which  X = 0 corresponds  to  the 
uniform  distribution,  R(v)(0,_£),  on  Sffl.  For  fixed  X = Xj  > 0,  the 
Neyman-Pearson  Lemma  gives  the  most  powerful  size  a test  as  rejecting 
the  null  hypothesis  for  large  values  of  the  likelihood  function 


n 

n r 
i=l 


(v)^1i'si)  = OV  exp(xi£'  l 

i=l 
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But  this  likelihood  function  is  a monotone  increasing  function  of 
= rn  _£'xq,  i = l,...,n.  Hence  the  uniformly  most  powerful 
(UMP)  size  a test  of  HQ:  X = 0 against  : X > 0 rejects  Hq  if 

rn  I’io  > k 

where  k is  chosen  so  that  P(rn  _£'xq  ^ k I ^ = ^)  = a. 

Equivalently,  rejecting  Hq  for 

VI)  = m n_l/2  h rnl'io  > c 

where  PCY^C^)  > c | X = 0)  = a is  the  UMP  size  a test  for  testing 
Hq:  X = 0 against  H^:X  > 0 when  _£  is  known. 

4.2.2  Uniformity  versus  Langevin  Distribution  for  Mean  Direction  Unknown 

If  _£  is  unknown,  we  will  show  that  rejecting  the  null  hypothesis  of 
uniformity,  Hq:  X = 0,  for  rn  sufficiently  large  provides  the  uniformly 
most  powerful  invariant  (UMPI)  size  a test  of  uniformity  against  the 
langevin  alternative,  H^:  X > 0. 

First  we  fix  X = X^  > 0.  Then  using  the  Neyman-Pearson  Lemma,  an 
invariant  test  under  orthogonal  transformations  is  found  by  rejecting 
uniformity  for  large  values  of 

L 11  r(v)< W dlj(v)(£> 

sm  i=l 

= Avn(X!)  J_  exp{X  5*  l u } U(v)(u_) 

Sm  i=l 

= AyD( Xt ) / exp{xl  *n£%}  dU(v)U) 

bm 
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(*1>  Av(Xrn)- 


Since  for  x > 0, 


t-  A (x) 
dx  v 


» 1 (x2/4)k 

= -d  l 

dx  k=0  (v+1).  k! 

k 

2 oo  1 (x2/4)k 

= — I 

X k-1  <v+l)k  Oc-1) ! 

x °°  1 (x2/4)k 

2 k=0  (vfl>k+1  k! 

x “ 1 (x2/ 4)k 

I 

2( v+1)  k=0  (v+2)k  k! 

x 

= A , , (x) 

v+1 

m 


> 0, 

we  see  that  Av~n(X1)  Av(X1r£1),  \l  > 0,  is  an  increasing  function  of 
rn.  Thus,  rejecting  uniformity  when  rn  > k where  P(rn  > k j X = 0)  = a 
provides  the  UMPI  size  a test  of  Hq:  X = 0 versus  the  Langevin 
alternatives  Hp  X > 0. 

In  section  3.3  we  showed  that  = m n 5^  rn»  Hence 

rejecting  Hq.'  X = 0 for  Y^(_Xq)  > c with  P(Y^(5q)  >c 
the  UMPI  size  a test  against  H^:  X > 0. 


X = 0)  = a gives 
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4.2.3  Uniformity  versus  the  Generalized  Cardioid  Distribution 

Here  we  derive  a locally  most  powerful  (LMP)  size  a test 
prescription  for  testing  uniformity  aginst  alternatives  of  the  form 

(j>p(_u'_£)  = 1 + ap  a£v)(u'_0  p = 1,2,3,... 

in  which  jj_  is  known.  In  order  to  insure  that  <j)p(u_'_£),  p = 1,2,3,...  is 
a density  on  Sm  we  require  0 < ap  < 1/A^V^(1).  Since  for  p = 1,  <f>^(u_'_£) 
is  the  cardioid  distribution  we  shall  call  <j>p(u_'_0 , p = 1,2,3,...,  the 
generalized  cardioid  distribution. 

A test  of  uniformity  versus  the  alternative  <|>  (u'jj)  may  be  written 

in  the  form  of  testing  Hq : ap  = 0 against  HQ:  ap  > 0,  p = 1,2,3 To 

derive  the  LMP  test,  we  shall  need  regularity  conditions  which  permit 
the  interchange  of  differentiation  and  integration.  Then  the  slope  of 
the  power  function  may  be  calculated  as 


I 3 

where  VCu^ , . . . | ap)  = ---  In  II  <|>  is  the  score  function  and  C 

p i=l 

is  the  critical  region  consisting  of  those  points  for  which 

uniformity  is  rejected. 

Verification  of  interchanging  differentiation  and  integration 

requires  n . ♦_(§>,)  to  be  bounded  for  0 < a < 1/A(v)(l), 

3a  i P 1 p p 

P 


P = 1,2,3 


Hence,  for  p = 1,2,3,..., 
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3a 


3 


P 


P 


< l a'^co  n {i  + y£v>(i» 


i=l  ^ j=l 


< n a£v)(1)  2n_l. 


Now  let  fQ  and  be  real-valued  integrable  functions  relative  to  a 
measure  p.  Then  an  extension  of  the  Neyman-Pearson  Lemma  (see  Fraser 
(1976))  gives  f^^)  > k fgCj)  as  the  critical  region,  C,  which  maximizes 
/ c f^i)  dy(^)  subject  to  fQ  £q(z)  duty)  = a. 

Applying  this  extended  version  of  the  Neyman-Pearson  Lemma  to  the 
slope  of  the  power  function,  we  find  the  LMP  size  a test  will  reject 
uniformity  for  large  values  of  the  score  function  evaluated  at  a^  = 0. 
That  is,  the  LMP  test  of  size  a has  critical  region 


c ~ { (u_i , . . . ..u^  : v(jLi » • • • »iinlap^  _Q  > kl 


P 


where 


>k)  = a,  p = 1,2,3 


The  score  function  evaluated  at  ap  = 0,  p = 1,2,3 


n 


ap=0  8a 


in  n $ pCi’Ui) 
i=l 


P 
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n 


A(v)(£'u  ) 

__2 Z-li- 

(v), 


i=l  l + a A'  "(E’u.) 

P P _1 


0 


n 


I 

i=l 


For  p = 1,2,3,...,  set  Tp  = n-l/2  l ± We  may  use  Tp 

directly  as  the  test  statistic  and  reject  uniformity  for  large  values  of 

T . Similar  to  the  results  for  the  mean,  variance  and  asymptotic  null 
P 

distribution  of  Yp(_£),  we  find  that  under  uniformity  Tp  - N(0,a£v^(1)) 
as  n -*■  °°.  For  a continuous  non-uniform,  rotationally 

symmetric,  unimodal  alternative,  Tp  - N(y,c>  ) as  n » where 


y = n^2fpA^v^(l)  and  a2 


l f2rl(v)Cp»P»2r) 
r=0 


Since  for  p = 1,3,5,..., 


n 


I 

i=l 


Vl>- 


we  may  use  the  results  already  derived  for  Yp(_£).  Thus  a LMP  size  a 

test  rejects  Hq  : ap  = 0,  p = 1,3,5 for  Yp(jj)  > k with  P(Yp(_£)  >k  | 

ap  = 0)  = a. 

As  a special  case,  we  observe  that  <j>|(_£'_u_)  = l + a^  A^v^(_£'_u), 

_£,u_  e Sm,  0 < at  < 1/m  is  the  cardioid  distribution  on  Sm.  Consequently 
for  given,  rejecting  Hq:  a ^ = 0 whenever  Y|(_£)  > k with  P(Yj(_£)  > k j 
a^  = 0)  = a provides  the  LMP  size  a test  against  the  cardioid 


alternatives  H^:  0 < a^  < l/m. 
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4.2.4  Uniformity  versus  other  Selected  Distributions 

With  only  slight  modifications,  the  derivation  in  subsection  4.2.3 
may  be  adapted  to  derive  a LMP  test  of  uniformity  against  the 
alternative 


+S(b(l'u.)  - 1 + b ( I ap  A^v)(i'i)t 

for  E,  given  and  0 < b < 1.  In  order  for  <f>s  b0  be  a density  on 

S , the  coefficients  a , p = 0,...,s,  must  be  sufficiently  small, 
m P 

The  testing  problem  may  be  formulated  as  testing  Hq:  b = 0 against 
Hp  0 < b < 1.  The  score  function  V (uj  , • • • | b)  evaluated  at  b = 0 

becomes 


(v) 


I b>  lb=o'4oa» 

Then  with  T = n-1/2  Y . a£v)U'u,),  a LMP  size  a test  may  reject 
p u l P 1 

uniformity  for  \ p apn^2  Tp  > k,  p = 0,...,s,  in  which  P(£  papn  ^p  > 
k | b = 0)  = a. 

If  the  even  numbered  coefficients  of  <(>s  ^(^'u)  are  zero,  define 


<f>s,b(l'a>  = 1 + M I a2p+l  "2p+l 

p=0 


A2d+1^'-^  1 * 


Hence, 

s 

V ( Hi  > • • • >iin  I I = C2p+1  Y2p+1^-^ 

where  c2p+1  = a2p+1  tt  n1/2  52J+1  , p = 0,l,...,s.  Then  the  critical 
region  for  a LMP  size  a test  of  Hq:  b = 0 against  Hji  0 ( b < 1 rejects 


for 
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I co  +1  Yo  +i (5)  > k 

p-o  2p+1  2p+1  ' 


such  that 


P(  jQ  c2p+l  > k I b ■ 0)  ' - 

In  addition,  Y^Xq)  shares  one  further  optimal  property  with  the 
Rayleigh  test.  Using  the  conventional  polar  angle  notation  as  in 
section  1.4,  we  may  project  the  bivariate  normal  distribution  N2(_n,E) 
onto  the  circle  S2  to  obtain  the  density  of  the  offset  normal 
distribution  on  S2  as 

gA(0)  = exp(-l/2  A2)  + Asin0  <f>(Acos0)  $(Asin0),  (dU(0)) 

0 < 0 < 2tt , where  $(•)  and  $(•)  are  the  standard  normal  density  and 

2 ”1 

distribution  functions  respectively  and  A = _n'£  _n.  Bhattacharyya  and 

Johnson  (1969)  have  shown  that  the  Rayleigh  test  and  hence  Y^(jcq)  is 
locally  most  powerful  invariant  for  testing  uniformity,  Hq:  A = 0, 
versus  the  offset  normal  alternatives,  H^:  A 4 0. 
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4.3  Power  Computations 

We  have  noted  optimal  properties  of  the  test  statistics  Yp(_§)  and 
Yp(xQ)  for  the  Langevin,  generalized  cardioid,  offset  normal  and  other 
selected  distributions.  For  the  alternatives  concerning  Saw's 
distribution.  Saw's  hyperspherical  t distribution,  the  Cauchy 
distribution  and  the  wrapped  normal  distribution,  tests  involving  Yp(_£) 
or  Yp(jjcQ)  for  uniformity  with  optimal  properties  are  presently  not 
known. 

However,  foreknown  the  statistic  Y^(_£)  appears  to  have  good  power 
for  the  group  of  alternatives  listed  above.  For  these  alternatives, 
there  seems  to  be  little,  if  any,  gain  in  power  by  considering  either  a 
linear  combination  of  Y^(_£)  and  Y^(_£)  or  a linear  combination  of  their 
squares. 

For  large  samples  and  _£  known,  we  may  reject  uniformity  when  Yp(jj) 

> c where  c = zq(  TjA[  ( l)  ) and  zq  is  the  upper  100a  percentage  point 
of  the  standard  normal  distribution.  For  large  samples  the  power  is 
then  given  by 


PHl(Yl(i-)  > c)  = l - $((c-Ul)/0l) 

where  <K  • ) is  the  standard  normal  distributin  function  while  and 
are  defined  in  Theorem  3.4.1.  In  the  following  tables  we  present  powers 
of  tests  rejecting  uniformity  for  Y^Cjj)  > c against  the  alternatives 
mentioned  above.  Rejecting  for  Y^(_£)  > c seems  to  perform  well  for  the 
Cauchy  and  wrapped  normal  distribution  but  less  well  for  Saw's 


distribution. 
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Table  4.1  Powers  of  the  Test  Y^(J^)  > c,  a = .05, 
for  Testing  Uniformity  against  the  Cauchy  Distribution 

with  Parameter  X 


n = 15 


X 

m 2 

3 

4 

5 

10 

15 

0.1 

.135 

.164 

.191 

.217 

.336 

.442 

0.2 

.287 

.378 

.461 

.535 

.794 

.917 

0.3 

.499 

.650 

.762 

.841 

.983 

.999 

0.4 

.724 

.871 

.944 

.976 

1.000 

1.000 

0.5 

.897 

.976 

.995 

.999 

1.000 

1.000 

0.6 

.980 

.999 

1.000 

1.000 

1.000 

1.000 

0.7 

.999 

1.000 

1.000 

1.000 

1.000 

1.000 

X 

m 2 

3 

n = 25 
4 

5 

10 

15 

0.1 

.173 

.217 

.258 

.298 

.474 

.615 

0.2 

.407 

.535 

.641 

.727 

.939 

.989 

0.3 

.691 

.841 

.922 

.963 

.999 

1.000 

0.4 

.902 

.976 

.995 

.999 

1.000 

1.000 

0.5 

.985 

.999 

1.000 

1.000 

1.000 

1.000 

0.6 

.999 

1.000 

1.000 

1.000 

1.000 

1.000 
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Table  4.2  Powers  of  the  Test  Y^(_£)  > c,  a = .05, 
for  Testing  Uniformity  against  Saw's  Distribution 
with  Parameter  X 


n = 15 


X 

m 2 

3 

4 

5 

10 

0.1 

.122 

.126 

.128 

.130 

.133 

0.2 

.242 

.256 

.264 

.269 

.279 

0.3 

.408 

.435 

.449 

.458 

.478 

0.4 

.596 

.631 

.648 

.661 

.684 

0.5 

.762 

.802 

.818 

.830 

.846 

0.6 

.893 

.917 

.927 

.933 

.942 

0.7 

.963 

.974 

.979 

.981 

.984 

0.8 

.991 

.995 

.996 

.996 

.997 

0.9 

.999 

.999 

.999 

1.000 

1.000 

n = 25 


X 

m 2 

3 

4 

5 

10 

0.1 

.153 

.159 

.162 

.164 

.169 

0.2 

.340 

.360 

.371 

.378 

.393 

0.3 

.580 

.613 

.630 

.640 

.661 

0.4 

.798 

.828 

.842 

.850 

.867 

0.5 

.932 

.948 

.955 

.959 

.966 

0.6 

.985 

.991 

.992 

.993 

.995 

0.7 

.998 

.999 

.999 

.999 

1.000 

15 

.134 

.283 

.485 

.692 

.852 

.945 

.985 

.997 

1.000 


l_5 

.171 

.398 

.669 

.872 

.968 

.995 

1.000 
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Table  4.3  Powers  of  the  Test  Y^(_§)  > c,  a - .05, 
for  Testing  Uniformity  against  Saw's  Hyperspherical  t-distribution 

with  Parameter  X 


n - 15 


n = 25 
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Table  4.4  Powers  of  the  Test  Y^(_£)  > c,  a = .05, 
for  Testing  Uniformity  against  the  Wrapped  Normal  Distribution 

with  Parameter  X 


n = 15 


X 

m 2 

3 

4 

5 

10 

15 

0.1 

.134 

.161 

.187 

.212 

.329 

.438 

0.2 

.283 

.373 

.458 

.539 

.848 

.981 

0.3 

.499 

.665 

.798 

.894 

1.000 

1.000 

0.4 

.742 

.915 

.986 

1.000 

1.000 

1.000 

0.5 

.928 

.997 

1.000 

1.000 

1.000 

1.000 

0.6 

.995 

1.000 

1.000 

1.000 

1.000 

1.000 

n = 25 

X 

m 2 

3 

4 

5 

10 

15 

0.1 

.172 

.215 

.255 

.294 

.473 

.624 

0.2 

.405 

.537 

.650 

.745 

.974 

1.000 

0.3 

.670 

.865 

.952 

.987 

1.000 

1.000 

0.4 

.921 

.992 

1.000 

1.000 

1.000 

1.000 

0.5 

.994 

1.000 

1.000 

1.000 

1.000 

1.000 

CHAPTER  5 

GRAPHICAL  DENSITY  ESTIMATION 


5.1  Introduction 


To  this  point  we  have  been  concerned  with  testing  uniformity  on  Sm 
against  certain  rotationally  symmetric,  unimodal  alternatives.  Besides 
conducting  tests  of  uniformity,  in  practice  it  would  be  desirable  to  be 
able  to  suggest  from  which  particular  distribution  we  are  sampling. 

Here  we  present  a graphical  curve  fitting  aid  for  distributions  on 
Sm,  similar  to  that  for  the  Pearson  frequency  curves,  i.e.  Types  I-VII, 
as  given  by  Pearson  and  Hartley  (1966).  As  an  aid  to  determining  the 

O O 

particular  Pearson  curve,  a graph  of  0 = P3  / plotted  against  g = 

/ ^2 , where  pr  is  the  rc  moment  about  the  mean.  In  our  case,  a 

distribution  on  Sm  may  be  determined  from  a graph  of  f^v^  plotted 
against  f where  v = m/2  - 1 and  fpV^  is  the  pC^  coefficient  in  the 
Gegenbauer  transform  of  the  density  f(X_£'u).  Since  most  practical 
applications  dictate  data  in  two  or  three  dimensions,  we  report  curve 
fitting  graphs  for  m = 2 and  m = 3 in  Figures  5.1  and  5.2  respectively. 


5.2  A Graphical  Curve  Fitting  Aid 
Suppose  that  f(X^__u^)  is  a density  on  Sm  with  Gegenbauer  transform 


v 

L 

p=l 


’bV)  A^Ui’u? 
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The  coefficients  of  the  Gegenbauer  transforms  for  many  distributions  on 
Sm  are  given  in  section  2.3.  At  the  end  of  this  section  Figures  5.1  and 
5.2  give  graphs  of  fjV^  plotted  against  f^v^  for  m = 2 and  m = 3 
respectively.  Given  the  random  sample  ... ^ from  Sm  we  may  estimate 


With  v = m/2  - l determined  by  context,  we  henceforth  discontinue  the 
superscript  (v)  notation  for  f£v^  and  fpv^  . 

a 

The  quantity  fp  is  unbiased  for  fp,  p > l,  as  is  seen  from 

EV  (nApV)(l)Tl  ^ EA^Cl'Uj) 

" (nA^v)(l))"1  l f a£v)(_£'u,)  I f A(vj)(^*u.)dU(v)(u..) 

j=l  Sm  J r=0  J J 

= KV)(1»-1  l / [A^V)(1'U j)]2  dU^>(u.) 

j—  l m 

■v 

A 

The  second  moment  of  fp,  p > 1,  is  found  to  be 

Efp  = (n  aW(1))-2  E { ^ A^>U'u.)}2 

- (nA^v)a)r2  { l e(a^)(1-u  ))2 

j=l  p J 

+ I I EApV)(i.,Jii)  ea(v)(_c'_u.)} 
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- (nA<»><l)r2  1 I f [A'v)(ru.j)]2 

j — 1 am 

x l frADV)(i.'iii)  dU(v)(u_.) 
r=0  J J 

+ IUL  Ajv)(i'Ui)  l [L  ^v)(rui)  «<"><«*)] 

i*j  Sm  r=0  Sm 

* [/  A^U'u  ) l f^d’u  )]} 

Sm  r=0  J 

- (nA^CD)-2  ( l I l2ll  [A<*W)]2 

j=l  r=0  sm  J 

x A^>d’Uj)  dU(v)(u.) 

+ I I [fp/_  ApV)d'iii)  AJU)d’iii)  dU(v>(u.)] 

x ^Vc  APV)(^J)  A£V)(lV  du(v)^j)]} 

- (n^U))'2  ( l ? f2rI(v)(p,p,2r) 

j=l  r=0 


+ 11  (V^O)]2) 

= n~l  { [ApV>  (l)  ]~2  l f2rI(v)(p,p,2r)  + (n-l)  f*} 
r=0  v 

where  I(v)(p,q,r)  = A<v)d'u)  aJv)(^'u)  a[v)(^u)  dU(v)(_u)  is  given 
by  formila  2.2.69.  Denoting  the  variance  of  by  y^  we  have, 


= n 


-'([A^CD]-2 


P 

I 

r=0 


:2rI(v)(p,p,2r)  - fj 


}• 
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Then  by  the  central  limit  theorem  we  find, 

A — 

f - f . 

— £ E-  5 N(0 , 1)  as  n -*■  ®. 

Y 

Suppose  the  random  sample  u^ ,...  is  given.  In  practice,  the 
modal  vector,  J,  is  usually  unknown  so  that  we  may  estimate  1 by  Xq,  the 
sample  mean  direction,  as  defined  in  section  1.5.  We  may  then  use 


■(.iWai)'1  l A<v,(io’ii> 

p j=! 


to  estimate  fp. 

For  f(X£/_u)  a density  on  Sm  we  may  construct  a graph  of  f x versus 
f2  for  various  values  of  X.  In  section  2.3  expressions  are  given  for  f x 
and  f2  for  the  Langevin,  Saw's,  Saw's  hyperspherical  t,  Cauchy,  wrapped 
normal,  cardioid  and  point  distributions.  Figures  5.1  and  5.2  present 
such  graphs  of  f L versus  f2  for  m = 2 and  m = 3 respectively.  By 
computing  f*  and  f2  and  observing  which  graphed  curve  is  closest 
to  the  point  (f*,f*),  one  may  decide  upon  a distribution  to  use  to  fit 
the  data. 
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Selected  Distributions  with  m=2 
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f, 


Figure  5.2  Graph  of  ^ versus  f for 
Selected  Distributions  with  m=3 


CHAPTER  6 
SUMMARY 


A brief  introduction  to  directional  data  was  presented  in  Chapter 
l,  section  l.l  and  there  it  was  argued  that  the  main  importance  of 
dealing  with  methods  applicable  to  general  m-spheres  is  to  unify  present 
methods  now  treated  separately  on  the  circle  and  sphere.  Some  areas 
where  directional  data  methods  find  application  were  indicated  in 
section  t.2. 

In  section  1.3  we  defined  uniform  measure  on  the  m-sphere,  Sm,  and 
densities  relative  to  this  uniform  measure.  Several  useful 
distributional  theorems  were  presented  in  section  1.4  along  with  a 
catalog  of  common  distributions,  namely,  the  Langevin  distribution. 

Saw's  distribution,  Saw's  hyperspherical  t-distribution,  the  circular 
cardioid  distribution,  the  circular  point  distribution  and  the  circular 
wrapped  normal  distribution.  Where  appropriate,  certain  properties  of 
these  distributions  were  recorded. 

Section  1.5  pointed  out  that  the  uniform  distribution  in 
directional  data  occupies  a more  prominent  position  than  does  its 
uniform  counterpart  in  the  linear  case.  Indeed,  under  mild  conditions, 
the  resultant  length,  r^,  and  the  mean  direction,  jcq,  are  independent  if 
if  only  if  we  are  sampling  from  the  uniform  density  on  Sm.  In  addition, 
the  uniform  distribution  appears  as  the  limiting  distribution  in  the 
central  limit  theorem  on  the  circle. 


105 


106 


Several  well  known  tests  applicable  to  Sffl  for  testing  uniformity 
were  stated  in  section  1.6.  These  tests  include  Rayleigh  s test,  a 
generalized  version  of  Agne's  test  and  a generalized  version  of  a class 
of  tests  given  by  Beran  (1968,  1969). 

Having  surveyed  the  current  tests  for  uniformity  in  section  1.6,  we 
presented  the  foundation  for  density  representation  on  Sm  in  section 
2.1.  It  was  shown  that  Gegenbauer  polynomials,  defined  in  section  2.2, 
provided  a means  to  represent  densities  f(_x'jj)  which  have  rotational 
symmetry  and  which  are  unimodal. 

The  definition  of  Gegenbauer  polynomials  given  in  section  2.2  used 
a standardization  and  weight  function  proposed  by  Saw  (1984)  . Here  it 
was  seen  that  this  particular  standardization  and  weight  function 
simplified  many  common  results  including  the  integrated  square  of  a 
Gegenbauer  polynomial.  The  remainder  of  section  2.2  was  devoted  to 
listing  many  standard  results  for  Gegenbauer  polynomials  in  terms  of 
this  standardization  and  weight  function. 

Section  2.3  defined  the  Gegenbauer  transform  of  a function  and 
listed  in  Table  2.5  the  coefficients  of  the  Gegenbauer  transform  of  some 
common  densities. 

Due  to  the  complicated  asymptotic  distribution  of  Ajne's  test 
statistic  generalized  to  m-dimensions , Chapter  3 proposed  a new  test 
statistic  for  testing  uniformity  on  Sm.  This  new  test  statistic  was 
shown  to  have  a normal  asymptotic  distribution.  But  before  giving  this 
new  test  statistic,  section  3.2  defined  X(u^)  to  be  a standardized 
version  of  N(u) , the  number  of  sample  points  in  some  m-dimensional 
hemisphere,  and  defined  <S(_u_'_£)  to  be  the  signum  function.  The  mean  and 
covariance  were  computed  separately  for  the  functions  X(_u)  and  6(u_'Jj) . 
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As  a final  preliminary  result,  the  coefficients  for  the  Gegenbauer 
transform  of  x(u/ w) , the  covariance  between  X(u)  and  X( w_ ) , was  given. 

Section  3.3  defined  the  new  test  statistic  Yp(_£)  and  gave  a 
computational  form  for  Yp(_g)  as  well.  Under  the  null  hypothesis  of 
uniformity,  the  mean  and  covariance  for  Y (jj)  were  derived  along  with 
the  fact  that  Yp(_£)  has  an  asymptotic  normal  distribution.  The  exact 
distribution  of  Y^(j^)  was  presented  for  m = 3.  The  exact  distribution 
of  Y^Xq)  was  also  given  under  the  hypothesis  of  uniformity. 

For  non-uniform  rotationally  symmetric  unimodal  alternatives, 
section  3.4  listed  the  mean  and  covariance  of  Yp(jj_)  and  again  found  that 
Yp(_£)  has  a normal  asymptotic  distribution. 

Chapter  4 discussed  testing  uniformity  with  Yp(jj) . Section  4.1 
showed  a relationship  between  the  statistics  (Yp(_£)}  and  the 
generalization  of  Ajne's  statistic  as  given  by  Saw  (1984). 

Optimal  properties  of  certain  tests  involving  the  statistics  Yp(_|) 
°r  Yp( Jq ) were  shown  in  section  4.2.  There  we  derived  UMP  and  UMPI 
tests  of  uniformity  versus  the  Langevin  distributin  for  mean  direction 
known  and  unknown  respectively.  We  gave  LMP  tests  against  the 
generalized  cardioid  alternative  and  other  selected  alternatives. 

For  some  alternatives,  tests  for  uniformity  involving  Yp(jj)  or 
V*o>  which  have  optimal  properties  were  not  known.  In  these  cases 
section  4.3  examined  the  power  of  Y^( %)  which  seemed  to  perform  well  or 
moderately  well  against  these  alternatives. 

In  the  course  of  these  investigations,  it  was  found  desirable  to  be 
able  to  suggest  from  wich  particular  distribution  a sample  was  taken. 
Thus  a graphical  curve  fitting  aid  was  presented  in  Chapter  5. 
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Of  future  interest  would  be  to  more  closely  examine  the  power  of 
certain  tests  involving  Yp(_£)  or  Yp(xQ>.  For  instance,  if  we  reject  for 
large  values  of  c^Y^(_£)  + c^Y^Cj^),  then  to  choose  c^  and  C3  to  maximize 
the  power  involves  minimizing  the  function 

0°13  f Vl  + C3P3, 

[z  } 

°13  ° 0al3 

where  p^  and  are  the  respective  means  of  Y^(_£)  and  Y^(_^)  while  q°13 
and  a^3  represent  the  covariance  between  Y^(_£)  and  Y3 (_£)  under  the  null 
and  alternative  hypotheses  respectively. 

It  would  also  be  of  interest  to  perform  some  simulation  studies  to 
compare  this  new  test  statistic  with  the  generalized  version  of  Ajne's 
statistic  and  to  check  the  rate  of  convergence  to  normality  of  the 
distribution  of  Yp(_£). 


APPENDIX:  INDEX  OF  NOTATIONS 


Def iniendum 

Em 

_x’X 

Sm 

x_ 


D<"> 

V 
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£ 
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Def iniens  Page  Introduced 

Euclidean  space  of  m dimensions  l 

inner  product:  _x'^_  = ^ x^y^  l 

surface  of  the  unit  m-sphere:  (x_:  x.  e Em,_x'_x  = l}  l 

column  vector  l 

surface  "cap":  [x:  x_  e Sm,  _x'_^  > u}  4 

collection  of  all  surface  "caps": 
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almost  everywhere  with  respect  to  a measure  y:  except 

on  a set  of  y— measure  zero  7 
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m 
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Def iniendum 
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xi 
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d 
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covariance  matrix  of  the  random  vector  x 11 
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Langevin  density  on  Sm  12 

Av(g) 
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modified  Bessel  function  of  order  v 12 
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1 q 
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Fisher  density  on  the  3-sphere  13 

R(u)(X,j) 
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density  for  Saw's  distribution  on  Sm  14 
Saw's  distribution  on  Sm  14 
density  for  Saw's  hyperspherical  t-distribution  15 
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Saw's  hyperspherical  t-distribution  15 
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generalization  of  Ajne's  statistic  I on  the  circle  23 
generalization  of  Ajne's  statistic  N on  the  circle  23 
Ajne's  test  statistic:  I = 2ira  A(tt)  23 
Ajne's  test  statistic:  N = B(tO  23 
null  hypothesis  23 
alternative  hypothesis  23 
limiting  distribution  of  the  Kolmogorov-Smirnov 

statistic  D(n)  24 
Kolmogorov-Smirnov  statistic  24 
m-dimensional  analogue  of  Ajne's  statistic  A(tt)  25 
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Laplacian  operator 
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Bessel  function  of  the  first  kind: 
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